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INTRODUCTION 

During  the  past  fifty  years  much  attention  has  been 
paid  to  the  physical  equilibrium  existing  between  gas  and 
liquid  phases.  Preliminary  investigations  involved 
relatively  simple  two  component  systems  but  as  the  need 
for  knowledge  of  the  behavior  of  more  complex  systems 
became  apparent  in  the  field  of  industry,  the  investigations 
were  extended  to  include  more  complex  systems  such  as  are 
discussed  in  this  thesis. 

The  present  work  is  concerned  with  a  relatively 
important  section  of  the  field  of  phase  equilibrium  studies 
—  viz  the  equilibrium  between  two  phase,  gas-liquid 
hydrocarbon  systems.  The  importance  of  this  work  arises 
from  the  necessity  for  knowledge  of  pressure-volume- 
temperature-composition  relations  as  applied  to  the 
petroleum  industry.  Until  quite  recently  accurate 
calculations  involved  in  oil  well  production  and 
petroleum  refining  were  limited  to  regions  of  relatively 
low  pressure  and  temperature  --  i.e.  regions  in  which  the 
ideal  laws  of  Henry  and  Raoul t  were  relatively  accurate. 
Recent  advances  in  petroleum  technology  have  necessitated 
a  knowledge  of  the  behavior  of  naturally-occurring 
hydrocarbon  systems  under  conditions  of  temperature  and 
pressure  far  beyond  the  region  in  which  these  ideal  laws 
are  applicable. 

Comprehensive  studies  of  phase  equilibria  over  wide 
ranges  of  temperature  and  pressure  in  two  phase  systems 
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of  this  type  may  be  considered  to  have  started  about  1930# 
Early  studies  involved  commonly  only  two  component  systems 
and  investigations  of  oomplex  nat urally-occurring  systems 
are  limited  very  largely  to  the  past  three  or  four  years. 
The  increasing-  interest  in  this  field  is  evidenced  by  the 
increasing  number  of  articles  appearing  in  current 
literature. 

The  purpose  of  such  investigations  is  primarily  to 
provide  the  necessary  data  for  accurate  calculations 
concerning  the  behaviour  of  naturally-occurring  systems 
in  the  production  and  processing  of  petroleum.  Theoretical 
interest  also  pertains  to  these  studies  in  so  far  as 
data  are  provided  for  evaluation  of  coefficients  in 
empirical  and  theoretical  equations  of  state  and  for 
confirmation  of  the  kinetic -molecular  theory  and  theories 
pertaining  to  intermolecular  forces.  These  data  are 
available  from  the  pressure-volume-temperature  relation¬ 
ships  of  the  gaseous  state,  determination  of  which  is  a 
necessary  preliminary  procedure  to  two  phase  equilibrium 
studies. 

The  purpose  of  the  current  work  is  of  a  two-fold 
nature.  First,  it  is  hoped  that  the  information  derived 
herein,  will  be  of  value  in  the  proper  operation  and 
exploitation  of  the  crude  oil  field  existing  in  Turner 
Valley.  At  the  present  time  investigations  are  being 
pursued  (2)  in  the  development  of  a  new  plan  of 
operation  of  wells  in  the  Turner  Valley  field  based  on 
removal  from  each  well  of  the  oil  and  gas  content  of  a 
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certain  maximum  volume  of  reservoir  spaoe  per  unit  area 
per  unit  time.  Proper  application  of  this  plan  depends 
almost  entirely  on  an  accurate  knowledge  of  the  phase 
densities  of  the  two  phases  in  the  producing  horizon. 

It  is  hoped  that  the  present  work  will  provide  information 
necessary  for  the  calculation  of  these  phase  densities 
and  also  of  the  effect  on  these  densities  of  decreasing 
reservoir  pressure  with  time. 

Second,  the  validity  of  certain  suggested  equations 
of  state  may  be  subjected  to  confirmation,  or  suitable 
new  expressions  may  be  derived.  The  derived  data  will 
serve  to  increase  the  already  large  fund  of  information 
concerning  fugacities  and  thermodynamic  properties  of 
gaseous  hydrocarbon  mixtures. 

This  work  is  a  continuation  of  a  series  of 
investigations  grouped  under  the  general  heading 
"Phase  Equilibria  in  Hydrocarbon  Systems".  The  project 
has  been  under  investigat ion  in  these  laboratories 
since  1931,  at  which  time  the  studies  were  limited  to 
simple  two -component  systems.  Although  such  studies  could 
have  been  continued  almost  without  limit,  it  was  found 
desirable  to  extend  the  scope  of  the  work  to  include 
complex  polycomponent  systems  more  nearly  related  to  those 
found  in  practice. 

Studies  on  complex  naturally-occurring  systems  were 
initiated  by  Gillies  and  Hugill  (6)  on  a  natural  gas- 
naphtha  sample  from  a  well  located  in  the  Turner  Valley 
gas  oap  area.  Upon  discovery  of  a  crude  oil  zone  in  the 
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Turner  Valley  field  it  was  felt  necessary  to  extend  the 
work  to  a  system  consisting  of  crude  oil-natural  gas. 
Preliminary  studies  on  a  system  of  this  type  using  Turner 
Valley  limestone  zone  orude  oil  and  Viking  natural  gas 
were  begun  by  Davis  and  Boyer  (4)  and  continued  by  Boyer 
and  Morrison  (1)  on  a  system  of  Turner  Valley  natural 
gas  and  crude  oil,  both  from  the  same  well.  Sufficient 
information  was  obtained  to  establish  the  method  used  and 
the  aoouracy  of  the  analytical  system  involved.  In 
addition  to  the  equilibrium  measurements  reported  by  these 
authors  (1,6),  pressure-volume-temperature  measurements 
on  natural  gas,  and  application  of  Lewis’  fugacity  rule 
to  hydrocarbon  mixtures  have  been  reported. 

The  present  work  is  a  direct  continuation  of  the 
work  of  Boyer  and  Morrison  with  a  view  to  investigating 
the  effect  on  the  two  phase  gas-liquid  equilibrium,  of 
variation  in  pressure,  temperature,  and  gas: oil  ratio. 

In  addition,  the  pressure-volume-temperature  relations 
of  the  natural  gas  have  been  determined  as  a  necessary 
preliminary  step  to  the  main  invest igation.  The  variation 
of  the  fugacity  of  the  gas  with  pressure  has  been 
determined  graphically  and  comparison  made  to  that 
calculated  using  Lewis’  fugacity  rule. 


. 


t 


. 

( ’< 


* 


, 

'  « t 


-  u 

.  .  +A 

- 


LITERATURE 


In  recent  years  the  number  of  papers  pertaining  to 
phase  equilibria  in  hydrocarbon  systems  has  become 
increasingly  large.  Comprehensive  reviews  of  the  literature 
concerning  solubilities  in  simple  systems  appear  in 
theses  by  Johnson  (8)  and  Piercey  (17).  Gillies  and 
Hugill  (6)  extended  this  bibliography  to  papers  concerning 
more  oomplex  hydrocarbon  systems.  Boyer  and  Morrison  (l) 
reviewed  the  literature  pertaining  to  pressure-volume- 
temperature  measurements  and  also  equilibrium  studies  of 
recent  date.  It  is  the  intention  in  this  review  to  refer 
only  to  those  papers  of  particular  significance  in  the 
present  work  and  to  those  papers  published  during  the  past 
year. 

Mention  might  be  made  of  the  work  of  Jessen  and 
Lightfoot  (7)  on  pressure -volume -temperature  data  for 
n-butane,  n-pentane  mixtures  and  for  the  pure  substances 
at  30°C.  Since  the  method  used,  involved  expansion  of 
a  known  volume  of  gas  into  a  larger  volume  and  measurement 
of  the  resultant  pressure,  it  is  felt  that  these  data  are 
not  as  reliable  as  those  of  Sage  et  al  referred  to  by 
Boyer  (1). 

Eielsen  Qf)  has  discussed  the  phenomenon  of  retrograde 
condensation  from  the  kinetic  molecular  viewpoint.  As 
pressure  is  increased,  the  Van  der  Waals  forces  of  the 
light  molecules  in  the  gas  phase  become  relatively  more 
important  than  similar  forces  tending  to  hold  the  heavier 
molecules  in  the  liquid  phase.  These  forces  between 
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heavier  molecules  in  the  liquid  phase  are  also  weakened 
by  solution  of  light  molecules  in  the  liquid  phase.  The 
result  is  transfer  of  heavy  moleoules  to  the  gas  phase 
with  increasing  pressure. 

Reference  has  already  been  made  (1)  to  the  work  of 
Kay  on  equilibrium  studies  in  the  systems  ethane -normal 
heptane  and  ethane-normal  butane.  This  work  has  been 
continued  (11)  to  include  the  normal  butane -normal  heptane 
system.  By  comparison  of  equilibrium  constants  for  butane 
in  ethane  and  heptane  and  for  heptane  in  ethane  and  butane 
the  author  concludes  that  the  critical  temperature  of  any 
binary  system  involving  two  of  these  components  is  closely 
proportional  to  the  sum  of  the  products  of  weight  fractions 
and  critical  temperatures  of  the  pure  constituents.  In 
all  three  binary  systems,  the  critical  pressure -compos it ion 
graph  shows  a  maximum  which  is  related  to  the  absolute 
molecular  weight  and  the  difference  in  molecular  weights 
of  the  constituents.  The  first  members  of  the  paraffin 
series  exert  a  greater  effect  than  do  higher  members. 

Sage,  Lacey  et  al  (3,  16,  20,  22)  have  continued 
their  extensive  studies  on  two-component  hydrocarbon 
systems,  and  subsequent  calculation  of  thermodynamic 
properties  of  these  systems.  Sage  and  Lacey  (19)  have 
attempted  to  correlate  equilibrium  constants  for  methane 
in  a  variety  of  hydrocarbon  systems  with  temperature, 
pressure,  molecular  weight  and  chemical  nature  of  the 
components  other  than  methane.  Ho  account  has  been  taken 
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of  the  variation  in  the  ratio  of  methane  to  less  volatile 
constituent.  The  results  are,  at  best,  qualitative  but 
provide  a  good  correlation  of  existing  data  in  a  variety 
of  systems.  The  authors  advise  against  unrestricted 
application  of  these  correlations  to  systems  not  strictly 
paraffinic  in  nature. 

Sage,  Hicks  and  Lacey  (21)  provide  a  correlation  of 
data  on  equilibrium  constants  for  lighter  paraffins  from 
methane  to  heptane  in  naturally-occurr ing  systems.  The 
treatment  of  methane  and  ethane  and  a  less  volatile 
constituent  is  fairly  detailed  and  the  results  are 
considered  reliable.  Data  for  all  other  hydrocarbons 
considered,  were  supplied  almost  entirely  from  the  work  of 
Katz  (9)  and  based  on  the  behavior  of  ideal  solutions  in 
low  pressure  regions.  Due  to  lack  of  available  data  the 
paper  does  not  discuss  the  effect  of  variation  in 
relative  amounts  of  the  different  constituents.  Discussion 
of  the  effect  of  constituents  lighter  than  pentane  on  the 
equilibrium  constants  for  lighter  hydrocarbons  is  omitted 
also. 

Shiah  (23)  has  made  use  of  a  simple  integration  of  the 
Clapeyron  equation  to  correlate  all  existing  plots  of 
equilibrium  constants  vs  temperature  and  pressure.  He 
shows  that  a  plot  of  the  natural  logarithm  of  the 
equilibrium  constant  vs  the  reciprocal  of  the  absolute 
temperature  is,  for  practical  purposes,  a  straight  line 
for  systems  on  which  data  are  available.  On  plotting 
reciprocals  of  temperature  for  different  substances  against 
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the  reciprocal  of  the  temperature  at  which  a  reference 
substance  has  the  same  equilibrium  constant,  a  straight 
line  is  obtained.  To  determine  the  equilibrium  constant 
at  any  temperature  for  any  hydro oarbon,  one  needs  only  to 
find  from  this  plot,  the  temperature  at  which  the  reference 
hydrocarbon  (pentane  for  convenience)  has  the  same 
equilibrium  constant.  From  a  plot  of  equilibrium  constant 
vs  temperature  for  the  reference  hydrocarbon  (pentane), 
the  value  of  the  equilibrium  constant  may  be  determined. 

The  graph  may  be  adapted  for  use  in  complex  systems  where 
the  average  molecular  weight  of  the  hydrocarbons  is  known. 

Dodge  and  Newton  (5)  develop  a  differential  equation 
relating  pressure  and  composition  of  both  liquid  and  vapor 
phases  in  equilibrium  at  constant  temperature.  Identical 
developments  are  made  using  free  energy  and  fugacity.  The 
equation  may  be  integrated  on  the  assumption  of  ideal 
conditions,  and  applications  to  specific  cases  are  shown. 
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THEORETICAL 


Original  considerations  of  two  phase  gas-liquid 
systems  at  equilibrium  were  based  on  the  ideal  laws  of 
Raoult  and  Henry*  Both  these  laws  assume  ideal  liquids 
and  gases,  i.e.  the  ohange  in  the  free  energy  (  &  F) 
is  proportional  to  the  logarithm  of  the  molar  concentrations, 
and  the  vapor  pressure  is  independent  of  the  total  pressure 
on  the  liquid.  At  low  pressures  they  apply  reasonably 
well  to  hydrocarbon  systems. 

That  liquid  and  gaseous  solutions  are  not  ideal  has 
been  long  known  and  many  devices  such  as  activity  and 
fugacity  have  been  developed  as  a  result  (12). 

That  vapor  pressure  is  a  function  of  total  pressure 
on  a  liquid,  was  shown  by  Poynting  (18).  He  showed  from 
thermodynamical  considerations  that  - 

(ggl)  r(l2) 

TO2'T;tX 

Where  P^  =  vapor  pressure  of  liquid 
Pp  *  total  pressure  on  liquid 
Tz  =  partial  molal  volume  of  liquid 
Yi  £  partial  molal  volume  of  vapor 

In  words,  the  increase  in  vapor  pressure  of.  a  liquid  with 

increase  in  total  pressure  on  the  liquid,  is  proportional 

to  the  partial  molal  volumes  of  liquid  and  vapor.  The 

equation  may  be  integrated  easily,  assuming  that  no 

solution  of  gas  in  liquid  occurs  with  increasing  pressure. 

Such  an  assumption  makes  the  integration  invalid  as  far 

as  natur ally-occurring  hydrocarbon  systems  are  concerned. 

For  application  to  the  present  case,  it  is  necessary 
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to  have  an  equation  for  each  phase  relating  composition 
and  pressure.  Johnson  (8)  develops  differential  equations 
of  this  type  for  two-component  systems  but,  before 
integration  can  be  carried  out,  it  is  necessary  to  have  an 
equation  of  state  for  each  phase.  Such  equations  of 
state  are  not  readily  available  at  the  present  time. 

Johnson  concludes  that  it  is  easier  to  measure  experimentally 
the  properties  of  the  individual  phases  and  set  up 
empirical  rules  than  to  evaluate  the  coefficients  of  an 
equation  of  state  for  each  phase.  In  the  present  studies 
on  poly-component  systems,  this  conclusion  appears  even 
more  applicable. 

In  representing  a  two-component  gas-liquid  system, 
the  customary  practice  is  to  make  use  of  G-ibb»s  Phase 
Rule  F=C-P*2.  By  restricting  one  of  the  variables, 
a  univariant  system  is  obtained  which  is  readily 
represented  graphically.  However  similar  representation 
of  a  polycomponent  system,  as  in  the  present  case,  is  more 
difficult.  In  view  of  the  difficulties  entailed  in 
developing  a  suitable  equation  of  state  for  each  phase, 
it  has  become  customary  to  represent  equilibrium  conditions 
by  the  so-called  "equilibrium  constants".  The  constant  is 


where  0-,^=  mole  fraction  of  component  1 
in  the  gas  phase. 

C-^  =  mole  fraction  of  component  1 
in  the  liquid  phase. 

A  true  constant  would  be  expected  if  both  phases  behaved 
ideally.  Since  this  is  not  the  case,  a  variation  in  K 
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with  pressure  at  oonstant  temperature  is  not  unexpected. 

Only  by  replacing  concentrations  by  activities  or 
fugacities,  can  a  constant  ratio  be  expected.  The 
constant  K,  as  defined,  would  be  expected  to  follow  the 
general  trend  of  the  constants  reported  by  Katz  and 
Hadkmuth  (9)  and  other  similar  investigations.  In  the 
case  of  complex  hydrocarbon  systems,  it  may  be  predicted 
that  K  for  methane  would  decrease  regularly  with 
increasing  pressure,  due  to  solubility  of  methane  in  the 
oil,  and  vaporization  of  the  oil,  until  a  value  of  1  is 
reached  at  the  critical  pressure  where  the  system  becomes 
single  phase.  Constants  of  less  volatile  components  should 
show  an  initial  decrease,  followed  by  an  increase  to  1, 
due  to  transfer  of  less-volatile  components  to  the  gas 
phase.  The  physical  result  of  this  effect  of  Poynting’s 
rule  gives  rise  to  the  so-called  "region  of  retrograde 
condensation". 

Attention  is  drawn  to  the  discussion  by  Katz  and 
Kurata  (10)  on  the  proper  application  of  the  term 
"retrograde  condensation".  In  this  region,  at  constant 
temperature,  an  increase  in  pressure  results  in  an 
increased  volatilization  of  liquid  into  the  gas  phase. 

It  is  in  this  region  that  the  equilibrium  constants,  K, 
for  the  liquid  components,  having  passed  through  a 
minimum  value,  are  increasing  with  increasing  pressure. 

The  variation  of  phase  densities  with  pressure  is  as 
would  be  predicted.  The  liquid  phase  density  shews  a 
continued  decrease  with  increasing  pressure  due  to  solubility 
of  lighter  components  from  the  gas  phase  and  transfer  of 
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heavier  components  to  the  gas  phase.  The  gas  phase 
density  shows  a  continued  inorease  due  to  the  normal  effect 
of  compression  of  the  gas,  transfer  of  heavier  components 
from  the  liquid  phase  and  preferential  solution  of  lighter 
components  in  the  liquid  phase.  At  the  critical  pressure, 
the  phase  densities  become  identical. 

While  it  is  always  found,  and  may  be  predicted,  that 
the  solubility  of  gas  in  liquid  increases  with  pressure 
on  the  system,  the  behavior  with  respect  to  changing 
temperature  is  not  so  simple.  For  ideal,  simple  systems, 
solubility  of  gas  in  liquid,  decreases  with  increasing 
temperature.  For  a  real  system,  this  is  not  necessarily 
true  and,  in  fact,  the  phenomenon  of  increasing  solubility 
with  increasing  temperature  might  be  expected  in  regions 
of  high  pressure  and  temperature.  The  result  of  this 
effect  on  phase  densities  is  a  rather  rapid  decrease  in 
liquid  phase  density  with  increasing  temperature.  Given 
the  correct  conditions,  a  simultaneous  increase  in  gas 
phase  density,  due  to  rapidly  increasing  volatility  of 
heavier  components,  may  occur. 

The  effect  is  common  and  has  been  demonstrated  in 
these  laboratories  by  Johnson  (8).  It  is  associated 
with  relatively  greater  transfer  of  the  normally  liquid 
components,  as  compared  to  transfer  of  gas  components,  to 
the  gas  phase.  The  effect  of  increasing  temperature  on 
the  position  of  the  region  of  retrograde  condensation  is 
to  reduce  the  pressure  at  which  the  phenomenon  occurs. 

Although  very  little  information  is  available  on  the 
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effect  of  variation  in  the  gas/oil  ratio  on  the  system 
at  equilibrium,  it  seems  reasonable  to  predict  that  an 


increased  gas/oil  ratio  would  tend  to  decrease  the  gas 


phase  density  due  to  dilution  of  the  heavier  components 
with  lighter;  simultaneously  the  liquid  phase  density  might 
be  expected  to  increase  due  to  more  complete  transfer  of 
more  volatile  components  to  the  gas  phase. 

As  has  been  mentioned  before,  the  determination  of 
pressure-volume-temperature  data  on  the  natural  gas 
sample  is  a  necessary  preliminary  step  to  the  equilibrium 
studies.  These  data  can  be  expressed  in  the  form  of 
empirical  or  theoretical  equations  of  state,  several  of 
which  have  been  proposed  elsewhere.  Any  such  equation  of 
state  is  limited  in  its  application  to  the  relatively 
narrow  limits  over  which  its  coefficients  are  valid.  A 
convenient  representation  of  pressure-volume-temperature 
data  lies  in  the  fugacity  (12).  Fugacities  were  determined 
on  the  Turner  Valley  natural  gas  sample  used  in  this 
investigation,  by  evaluation  of  the  expression  - 


7 

Where  f  =  fugacity 

a  s  residual  volume 
=  RT  -  V 

F~ 

and  P,  V,  T,  R,  have  their  usual  signif icance. 


The  residual  volume,  a,  is  defined  as  the  difference 
between  the  actual  volume  of  the  gas  at  pressure  P,  and 
the  volume  that  would  be  occupied  by  the  gas  if  ideal. 
The  equation  is  developed  from  the  fundamental  equation 
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dF 

3T 


Substitution  for  V  by  the  ideal  gas  law  leads  to  the 
simple  result 


dF  =  RT  d  InP. 


For  a  real  gas,  Lewis  (12)  suggested  the  use  of  the 
quantity  "fugacity" ,  f,  in  place  of  P  to  avoid  the  use 
of  a  complex  equation  of  state.  Thus  - 


dF  =  RT  d  Inf 


Or  AF  z  RT  lnfl 

r2 

As  P  — *  0 ,  the  gas  behaves  more  ideally  and  hence 


f  — >P  as  P— *0 


Combining  this  with  the  above  expression,  dF  - 

cTF  ”  7 


we  get 


RT  d  Inf  :  ?  dP 


Substituting  RT  -  a  for  V,  we  get 


Y~  p  p 

RT  lnf^  -  J  2  RT  dP  -  S 2adP 


Yo  Pi  P  Pi 

x  P 


Or  RT  In  f  =  RT  In  P  -  j'  adP 


o 


since  f— »P  as  P — * 


0 


P 

Thus  RT  In  f  z  “ f  adP 
r  o 


The  integral  was  evaluated  graphically  from  a  plot  of  a  vs  P 
and  hence  values  of  f  were  obtained. 

At  low  pressures,  it  is  possible  to  calculate  with 
reasonable  precision  the  partial  gas  pressures  in  a  gas 
mixture  of  known  composition  and  total  pressure  by 
Dalton’s  law.  However  in  regions  of  higher  pressures, 
Dalton’s  law  fails  to  apply  and  recourse  must  be  had  to 
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other  means.  Lewis  ( 12 )  has  proposed  a  useful  rule  of  the 

form : 

fx  =  %  f1°  and 

ft  =  *-  f2°N2  *- - 

Where  f^  z  fugacity  of  component  1  in 
a  gaseous  mixture. 

N^_  r  mole  fraction  of  component  1  in 
a  gaseous  mixture. 

f^°  z  fugacity  of  pure  component  1  at 

the  same  temperature  and  pressure 
^t  -  fugacity  of  the  mixture. 

Verification  of  this  rule  would  make  it  possible  to 

calculate  the  fugacity  of  a  gaseous  mixture  from  its 

composition,  knowing  the  fugacities  of  the  pure  components. 

The  fugacities  of  many  pure  substances  are  fairly 

accurately  established. 
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MATERIALS 

The  samples  of  natural  gas  and  crude  oil  used  in 
these  experiments  were  obtained  from  the  North  West 
Hudson* s  Bay  #6  well  in  Turner  Valley.  This  well  was  the 
farthest  south  producing  well  in  the  Turner  Valley  field 
and  is  looated  in  LSD  7,  Section  8,  Township  18,  Range  2, 
'West  of  the  3th  Meridian.  This  well  was  completed  on 
September  20,  1941  at  a  total  depth  of  823.5*.  Only  the 
upper  porous  zone  (7870*  -  7970*)  was  encountered.  In 
February,  1942,  when  the  samples  were  taken  the  well 
was  producing  at  the  rate  of  300  barrels  per  day  with  a 
gas/oil  ratio  of  1300  cu.ft.  per  barrel.  The  bottom-hole 
pressure  as  determined  on  October  19,  1941,  was  2377  psi 
and  the  bottom-hole  temperature  was  170°F. 

The  gas  sample  was  taken  through  a  fe1'  valve  on  the 
3"  fuel  line  from  the  separator.  The  oil  sample  was 
taken  through  the  bottom  valve  on  the  oil  level  gauge  on 
the  separator.  Eight  standard  gas  cylinders  and  three 
small  pressure  cylinders  were  used  to  collect  the  gas 
sample  and  a  fourth  small  cylinder  was  used  to  collect 
the  oil  sample.  The  oil  container  was  first  filled  with 
water  which  was  then  slowly  displaced  by  oil  from  the 
separator.  The  oil  sample  was  taken  continuously  at  the 
separator  pressure  at  such  a  rate  that  the  cylinder  was 
completely  filled  just  as  the  last  gas  cylinder  was 
filled.  The  gas  cylinders  were  flushed  three  times  by 
building  up  the  pressure  in  them  to  200  psi  and  then 
blowing  off  to  atmosphere.  Gas  samples  were  taken  every 
half  hour.  During  the  period  of  sampling,  the  well  was 
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produoed  at  a  steady  rate.  Sampling  was  done  by  the 
field  staff  of  the  Petroleum  and  Natural  Gas  Conservation 
Board,  supervised  by  Harrison  H.  Heisler. 

On  receipt  of  the  samples,  they  were  purified  in  a 
similar  manner  to  that  described  by  Boyer(l).  The  gas 
sample  was  scrubbed  to  remove  hydrogen  sulfide  and 
mercaptans  by  bubbling  through  40%  KOH  in  a  glass  column 
packed  with  short  lengths  of  glass  tubing.  The  clean 
gas  was  stored  in  a  gas  holder  from  which  it  was  compressed 
through  a  Calcium  chloride  drier  into  the  low  pressure 
storage  system.  The  gas  liberated  from  solution  in  the  oil 
was  added  to  that  from  the  gas  cylinders  and  the  whole  gas 
sample  was  mixed  until  a  constant  analysis  was  obtained. 

After  the  gas  in  solution  in  the  oil  had  been 
liberated,  the  oil  was  displaced  by  water  and,  after  being 
scrubbed  in  the  same  caustic  column,  was  collected  and 
stored  over  water. 

The  gas  composition  was  determined  in  a  low 
temperature  fractional  distillation  apparatus  (6).  An 
average  of  four  analyses  on  the  original  gas  gave  a 
composition  of: 

Component  Volume  % 


Methane 

78.95 

Ethane 

10.90 

Propane 

4.85 

Butane 

2.19 

Pentanes  *■ 

0.60 

Nitrogen 

2.55 

The  density  of  the  gas  was  determined  by  weighing 
a  calibrated  bulb  containing  the  gas  at  known  temperature 
and  pressure.  The  average  of  18  densities  was  8.151x10  ^g/oc. 
at  25°C  and  76  cm.  The  average  deviation  from  this  value 
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was  1*5  parts  per  1000. 

Density  of  the  crude  oil  was  determined  by  weighing  a 
known  volume  of  oil  at  20°C  and  25°C. 

d20°  «  o.815>1  g/cc. 

4° 

d2^  z  0.8111  g/cc. 

4° 

The  mean  Molecular  weight  of  the  crude  oil  was  determined 
by  the  cryoscopio  method  (6,17)  and  found  to  be  174. j?  g. 

The  crude  oil  was  fractionally  distilled  and 
refractive  indices  of  successive  fractions  plotted  against 
their  densities  and  mean  molecular  weights.  The  results 
are  shown  in  Table  1  and  plotted  in  Figure  1. 
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Table  1 


Molecular  weight 

-  refractive 

index  -  density 

measurements 

on  crude  oil 

fractions. 

%  over 

Molecular 

weight  Refractive  index 

Density 

5 

I.3862 

0.6724 

6 

83.O 

1.3784 

9 

83.4 

I.3848 

10 

1.3900 

o. 6867 

12 

88.8 

1.3509 

12*6 

101.9 

1.3902 

14 

80. 

1.5559 

15 

1.3941 

0.6929 

15 

93*8 

1.3954 

16.8 

91.9 

1.3966 

18 

91.7 

1.4008 

15.2 

97 

1.5999 

20 

1.3981 

0.7050 

20.2 

98.3 

1.4021 

21 

100.6 

1.4024 

21.6 

96.1 

1.4038 

25 

100.6 

1.4060 

25 

101.0 

1.4074 

Figure  1 

Variation  of  Refractive  Index 
with  Density  and  Molecular  Weight 
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APPARATUS 

The  experimental  apparatus  is  identical  with  that 
described  in  detail  by  Boyer  (1).  It  consists  of  three 
main  parts: 

1.  High  pressure  storage  systems  and  equilibrium 
cylinder. 

2.  Low  pressure  analytical  system. 

3.  Crude  oil  storage  and  dispensing  system. 

The  high  pressure  system  consists  essentially  of  a 
drier,  a  gas  measuring  pipette,  a  deadweight  piston  gauge 
and  the  equilibrium  cylinder  and  pipettes.  Gas  from  the 
storage  tanks  is  admitted  to  the  high  pressure  system: 
by  means  of  a  compressor  and  passes  through  a  KOH  drier 
surrounded  by  an  air  bath  at  60°C,  into  the  measuring 
pipette  which  is  kept  at  60°C~  O.IJ?0  by  means  of  a 
thermostatically  controlled  oil  bath.  The  approximate 
pressure  in  the  pipette  or  in  the  drier  can  be  read  on  a 
Bourdon  gauge.  The  pipette  can  be  connected  to  the 
dead. ;eight  piston  gauge  and  to  the  equilibrium  cylinder. 

To  prevent  solution  of  the  gas  in  the  castor  oil  of  the 
piston  gauge,  a  mercury  U  tube  is  inserted  in  the  line. 

Level  of  the  mercury  in  the  U  is  shown  by  deflection  of  a 
galvanometer  needle  placed  in  a  circuit  which  terminates 
in  a  platinum  point  at  the  mercury  surface  and  in  the 
wall  of  the  U  tube.  Pressure  in  the  pipette  may  be  increased 
by  a  hydraulic  mercury  pump. 

The  equilibrium  cylinder,  a  complete  description  of 
which  is  given  by  Boyer  (1)  is  totally  immersed  in  a 
thermostatically  controlled  oil  bath  heated  by  two  lj?00 
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watt  heaters,  one  of  which  is  connected  to  a  relay  cirouit. 
The  cylinder  may  be  raised  from  the  bath  for  charging  and 
for  talcing  equilibrium  samples. 

To  prevent  possible  condensation  of  higher  hydrocarbons 
under  high  pressure,  all  high  pressure  lines,  including 
the  compressor,  drier  and  U  tube,  are  maintained  at  60°C 
by  suitable  heaters. 

The  equilibrium  cylinder  consists  essentially  of  a 
hollow  steel  cylinder  3  1/16"  I.D.  Copper  connections 
join  the  cylinder  to  gas  and  liquid  sample  pipettes  which 
may  be  shut  off  by  valves  located  at  each  end  of  the 
pipettes.  A  pump  assembly,  driven  by  a  drive  shaft  passing 
through  one  end  of  the  cylinder,  circulates  gas  through  the 
gas  sample  pipette,  which  is  located  in  a  vertical  position 
above  the  cylinder,  and,  simultaneously,  liquid  through 
the  liquid  sample  pipette  which  is  located  in  a  vertical 
position  below  the  cylinder.  For  this  series  of 
experiments,  a  steel  half  cylinder  3"  long  and  2  13/16" 
in  diameter  was  fitted  into  the  equilibrium  cylinder 
between  the  pumps  and  the  opposite  end  plate.  This  plug 
reduced  the  volume  of  materials  needed  for  an  experiment 
and  permitted  the  use  of  much  lower  gas/oil  ratios  than 
•would  otherwise  have  been  the  case. 

The  analytical  system  is  unchanged  from  that 
described  by  Boyer  (1).  It  consists  essentially  of  a 
condenser  train,  a  gas  storage  and  measuring  system,  and 
the  low  temperature  fractional  distillation  column.  The 
apparatus  can  be  completely  evacuated  by  means  of  a 
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two-stage  mercury  diffusion  pump  backed  by  a  Iiyvao  oil 
pump.  The  sample  pipettes  can  be  connected  to  a  standard 
steel  connector.  The  first  condenser  is  removable  and  is 
connected  to  the  apparatus  by  a  ground  glass  joint.  The 
second  condenser  is  surrounded  by  a  mixture  of  dry  ice  and 
aoetone  and  the  third  condenser  is  cooled  by  dry  ice- 
acetone  under  vacuum.  These  mixtures  maintain  the  condensers 
at  -79°C  and  -112°C  respectively,  which  temperatures  are 
read  by  means  of  double  junction  copper-constantan 
thermocouples.  A  small  graduated  tube  of  about  5  cc.  capacity 
is  joined  to  the  condenser  train  by  means  of  a  thin-walled 
capillary.  This  tube  is  for  trapping  the  liquid  retained 
in  the  condensers  when  the  gas  phase  sample  is  expanded. 

The  liquid  is  distilled  into  this  tube  at  the  temperature 
of  liquid  air  and  the  capillary  is  sealed  with  an 
oxygen-gas  flame.  A  third  condenser  in  the  train  can  be 
surrounded  by  liquid  air  at  the  proper  time. 

The  condenser  train  is  connected  to  the  gas  measuring 
system  by  means  of  a  manually  operated  Toepler  pump.  The 
measuring  system  consists  of  six  1  litre  flasks  inside 
an  air  bath.  Any  or  all  of  these  flasks  can  be  connected 
to  the  Toepler  pump  and  to  a  constant  volume  manometer. 

The  low  temperature  fractionation  column  is  unchanged 
from  that  described  by  Gillies  and  Hugill  (6). 

The  crude  oil  is  dispensed  from  the  carboy  in  which 
it  is  stored  over  water  by  displacement  with  water  from  a 
res^voir.  It  is  dried  over  Calcium  chloride  and  measured 
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into  a  one  litre  flask.  After  weighing,  it  oun  he 
dispensed  into  the  equilibrium  cylinder  through  a 
capillary  tube  equipped  with  a  ground  glass  joint  to  fit  a 
standard  taper  male  connection  on  the  cylinder.  All  glass 
connections  are  lubricated  by  a  gel  of  mannitol  and 
glyoerol  in  which  the  oil  is  insoluble. 
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PROCEDURE 

Pressure -volume-temperature  relationships 

The  gas  sample  pipette,  evaouated  and  suspended  in  the 
oil  bath  at  the  proper  temperature,  was  filled  to  the 
desired  pressure  by  means  of  the  compressor  and  the 
hydraulic  mercury  pump.  Fine  adjustment  of  the  pressure 
was  made  on  the  deadweight  gauge  and  then  the  valve  on  the 
pipette  was  closed.  The  other  end  of  the  pipette  was 
connected  by  a  high  pressure  line  to  the  glass  analytical 
system.  With  all  necessary  parts  of  the  apparatus 
evaouated,  the  gas  sample  was  expanded  into  the  analytical 
system  and  transferred  into  the  bulbs  by  means  of  the 
Toepler  pump  until  a  residual  pressure  of  less  than  1  mm. 
of  mercury  could  be  read  by  means  of  a  McLeod  gauge.  The 
temperature  of  the  bulbs  was  brought  to  a  constant  value 
by  forced  circulation  of  air  and  the  pressure  in  the  bulbs 
measured  on  the  manometer.  Density  measurements  were 
made  on  this  gas.  Calculation  of  the  pres sure -volume - 
temperature  data  was  according  to  the  following  procedure: 

The  pres sure -volume  product,  PV,  of  the  compressed 
gas  at  temperature  T  is  represented  by 


(PV)  =  (pdw  *  PA)  ttgp 
14.696 


cc.  atmos. 


Where  P-^  z  Pressure  reading  on  deadweight  gauge  in  psi. 
P  =  Atmospheric  pressure  in  psi. 

VnT,  r  Volume  of  gas  pipette  in  cc. 

C-Jr 

The  pres sure -volume  product,  (P^V^^  of  the  expanded 
gas  at  23°C  and  76  cm.  of  mercury  is  expressed  by  - 
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(hh)1  '  fm 

76 


Va  + 


(V„  V,P)  P 


GP'  McL 


■m 


298.1 

T 


co.  atmos. 


Where  (P,  V,  z  PV  produo t  using  the  ideal  gas  laws  to 
1  oorrect  to  2^°  and  76  cm. 

Pm  z  pressure  of  manometer  in  cm.  of  mercury 

Va  r  volume  of  storage  system  in  co. 

Vr  z  volume  of  condenser  train  in  cc. 

VqP  =  volume  of  gas  pipette  in  cc. 

PmcL^  residual  pressure  on  McLeod  gauge  in  cm. 
of  mercury. 

T  z  absolute  temperature  of  storage  system. 
This  value  of  (PiV1)'1'  must  now  be  corrected  for 
deviation  of  the  gas  from  ideality  at  1  atmosphere. 
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r  1 


z  1 
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Where  (P,V,)  z  corrected  PV  product  at  29°C  and 
1  76  cm.  of  mercury. 

\  Z  deviation  factor  for  the  gas  at  1 
atmosphere. 


b 

AP 


r  deviation  factor  per  cm. 

-  76  -  P 

-  1  m 

is  then  plotted  against  P  (Figure  2) 


Equilibrium  Measurements 

The  equilibrium  cylinder,  previously  tested  for 
leaks,  was  set  in  place  in  the  oil  bath  and  valves  on 
both  pipettes  opened.  The  capillary  connection,  used  for 
charging  the  oil,  was  put  in  place  and  the  cylinder  was 
completely  evacuated. 


25. 

A  sample  of  oil  was  drawn  from  the  oil  storage  and 
chilled  in  ice.  After  weighing  the  flask  and  the  oil, 
the  latter  was  oharged  into  the  oylinder  through  the 
capillary  glass  connection  hy  allowing  atmospherio  pressure 
to  force  the  oil  into  the  evacuated  cylinder.  When  nearly 
all  the  oil  had  been  oharged,  the  vacuum  remaining  in  the 
oylinder  was  noted  on  a  manometer  attached  to  the  capillary 
glass  connection  and  then  natural  gas  was  admitted  to  bring 
the  pressure  up  to  atmospheric.  The  liquid  sampler  was 
removed  and  the  standard  steel  male  connection  on  the 
cylinder  was  quickly  closed  by  screwing  in  a  plug.  The 
weight  of  the  empty  flask  determined  the  weight  of  oil 
charged. 

The  volume  of  gas  at  25°C  and  76  om.  of  mercury  to 
give  a  desired  gas/oil  ratio,  was  then  calculated  and 
converted,  using  the  PV/P^V^  *  P  curve,  to  volume  at  60°C 
and  1^00  psi.  All  mercury  was  drained  from  the  gas 
measuring  pipette  and  the  gas  pressure  therein  was  raised 
to  exactly  1500  psi  at  60°C  as  measured  on  the  deadweight 
gauge.  This  gas  was  charged  into  the  cylinder  and  the 
pressure  in  the  pipette  restored  to  exactly  1500  psi  at 
60°C  by  draining  off  mercury.  The  actual  volume  of  gas 
added  at  1500  psi  and  60°C  was  found  by  weighing  the 
mercury  used  to  displace  it  from  the  pipette.  This 
procedure  was  repeated  until  the  desired  gas/oil  ratio  in 
the  cylinder  was  attained.  The  equilibrium  cylinder  was 
then  connected  directly  to  the  deadweight  gauge  and  the 
circulating  pumps  started.  When  a  constant  pressure  had 
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bean  reached  in  the  oylindar,  a  condition  of  equilibrium 
was  assumed  to  exist  and  the  pumps  were  shut  off. 
Equilibrium  was  usually  attained  within  an  hour  of  adding 
the  final  charge  of  gas,  particularly  if  the  pumps  had  been 
operating  during  the  gas  charging  operations.  The  cylinder 
was  rapidly  raised  from  the  bath  and  both  pipettes  closed 
off  at  once.  Pressure  was  released  from  the  cylinder  and 
the  pipettes  disconnected. 

Expansion  of  Gas  Phase  Sample. 

With  the  glass  apparatus  evacuated  and  the  proper 
temperatures  in  the  dry  ice-acetone  mixtures  around  the 
condensers,  the  phase  sample  was  expanded  slowly  into  the 
bulbs.  The  uncondensed  gas  was  transferred  to  the  bulbs 
to  a  residual  pressure  of  0.1  mm.  of  mercury  and  then  the 
third  condenser  was  surrounded  by  liquid  air.  This 
completely  removed  butanes  and  lighter  from  the  condenser 
train,  with  a  minimum  of  pentanes.  The  condensate  from 
the  liquid  air  trap  was  added  to  the  gas  in  the  bulbs  and, 
after  constant  temperature  had  been  attained,  the  pressure 
was  read  on  the  const ant -volume  manometer.  The  gas  was 
then  circulated  automatically  (I3)  until  it  was  thoroughly 
mixed.  Six  hours  was  considered  sufficient.  Duplicate 
density  measurements  were  made  and  samples  taken  for  low 
temperature  analyses. 

The  liquid  separated  in  the  dry  ice-acetone  condensers 
was  distilled  over  into  the  graduated  dilatometer  using 
liquid  air,  and  the  dilatometer  was  sealed  off.  Its 
weight  was  determined  and  the  volume  and  refractive  index 
of  liquid  at  20°C  measured. 
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Expansion  of  Liquid  Phase  Sample 

This  was  carried  out  in  identical  fashion  to  the  above 
To  prevent  violent  ebullition  of  gas  from  solution  in 
the  oil  at  very  low  pressures,  the  oil  was  stirred  by  a 
magnetic  stirrer  operated  manually  by  a  solenoid.  When 
the  last  of  the  butanes  and  lighter  components  had  been 
removed  to  the  bulbs,  the  condensers  were  lowered  and  a 
dry  ice-acetone  mixture  placed  in  the  cup  which  was 
soldered  to  the  pipette.  Thus  liquid  was  distilled  into 
the  pipette  and  it  was  thereby  washed  of  any  oily  residue. 
This  procedure  was  repeated  three  or  four  times.  Liquid 
air  was  then  placed  around  the  oil  condenser  and  all 
fractions  distilled  into  it.  The  condenser  was  removed 
and  fitted  with  a  ground  glass  cap.  The  weight  of  oil, 
its  density  at  20°C,  and  mean  molecular  weight  were 
determined.  An  A.S.T.M.  distillation  was  done  on  the  oil 
and  refractive  indices  determined  for  successive  fractions 
using  an  Abbe  refractometer. 


. 

' 

»  8 

- 

. 

* 

» 

,  ■  J ;  (tJL/01*^  .}  <1cr.lv;' 

.. 


I  i 


28 


RESULTS  AND  DISCUSSION 

The  results  are  presented  In  two  distinot  sections: 

1.  Pressure-volume-temperature  relationships  of  natural 
gas  from  N.W.  Hudson* s  Bay  #6,  and  comparison  with  other 
gases  of  a  similar  nature. 

2.  Two  phase  gas-liquid  equilibrium  measurements  on  the 
natural  gas-crude  oil  system. 

Pressure-Volume-Temperature  Relationships 
A  oomplete  series  of  PV  vs  P  values  was  determined 
for  this  gas,  at  60°C  in  the  pressure  range  1  -  190 
atmospheres.  One  point  was  determined  at  25  °C  and  two 
points  at  77°C.  The  experimental  results  are  presented 
in  Table  2.  Smoothed  values  at  60°C  obtained  from  a 
large-scale  plot  are  presented  in  Table  3. 

All  values  listed  in  Table  3  are  for  (PV)p^/(P]_V]_)2^°. 
To  convert  these  ratios  into  specific  values  of  PV/g  of 
gas  it  is  necessary  to  know  either  the  density  or  the 
molecular  weight  of  the  gas  for  a  particular  condition. 

Both  of  these  methods  were  used. 

The  average  of  18  densities  at  approximately  60  cm. 
of  mercury  and  25  °C  was  6.435  x  10  4g/cc.  The  average 
deviation  was  0.1 5%.  This  gives  a  value  for  (PV)2^o^0  cm. 

of  1226.84  cc,atm./g.  This  value  could  be  used  to 
calculate  the  molecular  weight  of  the  gas  but  inaccuracy 
due  to  deviation  of  the  gas  from  ideality  must  first  be 
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Table  2 


Experimental  values  of  PV/(PiV-|)  as  a  function  of 

Pressure  at  25°C,  60°C,  and  77°C  for 

fl.W.  Hudson* s  Bay  #6  Natural  Gas. 


P.  Atm. 

25°C 

60°C 

77  °c 

17.94 

1.0779 

1.0787 

34.95 

1.0441 

1.0444 

51.95 

68.97 

1.0137 

1.0129 

1.0108 

1.0105 

0.9736 

0.9737 

85.98 

0.9507 

0.9505 

102.99 

0.7274 

0.9299 

1.0147 

0.7304 

0.9303 

1.0129 

120.01 

0.9142 

0.9146 

137.02 

154.03 

0.9020 

0.9026 

0.8977 

0.8979 

0.8977 

0.9842 

171.04 

0.8985 

0.9832 

188.06 

0.9048 

0.9050 

Table  5 


Smoothed 

Values  of 

PV/(PnV-|)  from 

Curves 

through  Experimental 

Points  at  60°C. 

>.  atm. 

N.W.H.B.#6  A.O.I.  #2 

Royal it e  27 

0 

1.120^ 

1.1211 

1.1204 

20 

1.0776 

1.0789 

1.0811 

40 

1.0348 

1.0379 

1.0427 

60 

0.7943 

0.9991 

1.0078 

80 

0.9592 

0.9660 

0.9771 

100 

0.9532 

0.93?3 

0.9513 

120 

0.913? 

0.9175 

0.9317 

140 

0.9009 

O.9033 

0.9210 

160 

0.8972 

0.8954 

0.9168 

180 

0.9007 

0.8953 

0.9191 

200 

0.9052 

0.9279 

2?. 

There  are  several  ways  to  determine  the  deviation  of 
a  gas  in  the  regions  of  pressure  below  one  atmosphere.  The 
first  is  to  expand  a  known  volume  of  gas  into  successively 
larger  volumes  and  to  note  the  pressure  changes  at  constant 
temperature.  This  procedure  was  carried  out  by  Boyer  and 
Morrison  (1)  on  a  natural  gas  of  similar  composition  to 
the  one  used  in  the  present  experiments.  They  report  a 
linear  variation  of  PV  with  P  but  it  was  impossible  to 
determine  the  slope  of  the  graph  with  sufficient  accuracy 
for  the  purpose  in  mind. 

A  method  which  gives  considerably  greater  accuracy, 
is  to  extrapolate  the  graph  of  (PV)2^/ (^1^1)23  vs  p  from 
40  atmospheres  to  0  atmospheres.  In  this  range,  the 
graph  is  a  straight  line  and  extrapolation  is  very 
accurate  in  regions  of  low  pressure.  The  value  of  (PV)2^ 
at  60  cm.  of  mercury  is  accurately  known  from  density 
measurements,  and  from  the  slope  of  the  graph  of 
(PV)25/(PlVi)25  through  this  point,  it  is  possible  to  get 
an  absolute  value  of  (P]V.)oc.  Extrapolation  of  the  line 
to  zero  pressure  enables  one  to  calculate  an  absolute 
value  of  (P0V0)25  from  which  a  more  accurate  calculation 
of  the  molecular  weight  can  be  made  than  was  possible 
from  the  density  value  at  60  cm.  pressure. 

In  the  present  instance,  it  was  necessary  to  assume 
that  the  deviation  of  the  gas  from  ideality  in  the  region 
of  0  -  1  atmosphere  at  25°C  was  the  same  as  that  of  the 
Alberta  Oil  Incomes  #2  gas  used  by  Boyer  and  Morrison. 

This  gives  a  deviation  per  cm.  of  0.0000381,5.  This  was 
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used  to  oorrect  the  (P^V],)1  values  from  the  measured 

pressure  in  the  vicinity  of  60  om.  and  25°C  to  76  cm.  and 

25°C.  This  gave  the  corrected  value  of  (P,  V,  )oc.  From 

1  1 

the  slope  of  the  plot  of  (PV^q/CP^)  vs  P  in  the 
region  0  -  34  atmospheres  the  value  of  (P0V0)^q/(P1V1)2^ 
was  found  to  be  1.12047.  Since  Boyle*s  and  Charles*  Laws 
apply  exactly  to  ideal  gases,  they  can  be  applied 
accurately  to  real  gases  at  zero  pressure.  Thus  a  value 
of  (P0V0)2^/(P  1^)2^  is  readily  obtained  from  the  value  of 

(povo^60^plv1^23  usin&  the  factor  298.16/333 .16. 

This  yields  a  value  of  (P0V0)2^/(P1V1)2^  of  1.00276  as 
compared  with  1.00290  for  Alberta  Oil  Incomes  It  was 

considered  more  accurate  to  assume  a  correction  factor  for 
(PlV^)p  over  a  small  range  of  pressure  (16  cm.  approximately) 
and  from  this  corrected  value  to  calculate  the  slope  of 
the  (PV)2^/(P]_V]>)2^  vs  P  curve,  than  to  assume  the  same 
slope  as  existed  for  Alberta  Oil  Incomes  ^2  gas. 

It  is  obvious  that  the  value  of  (PV^cj/CP^V^)^  at 
one  atmosphere  is  1.0000.  From  this  value,  and  the  value 
calculated  for  the  slope  at  zero  atmospheres,  it  is 
possible  to  calculate  the  ratio  at  60  cm.  pressure.  Thus 

(^60v6o)25/(I>iVi)25  r  1  *■  16/76(0.00276)  =  I.OOO58. 

From  this  ratio  and  the  value  of  (P^qV^q)^  obtained  from 
densities,  we  get  (P^Vi)2^  =  1226.13  cc.  atm./g. 

Extrapolation  of  this  line  to  zero  pressure,  as  described 
above,  yields  a  value  of  (PqVo)2^  =  1229*51  cc. atm./g. 

This  value  yields  a  molecular  weight  of  19*91  g*  for  the 
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The  molecular  weight  of  the  gas  oan  also  he  found 
from  its  composition  as  determined  by  the  low  temperature 
fractionation  column.  The  result  by  this  method,  is  a 
moleoular  weight  of  20.4?  g.  This  is  a  considerably 
larger  error  than  was  found  by  Boyer  and  Morrison.  The 
discrepancy  can  tentatively  be  ascribed  to  inaocuracy  in 
experimental  technique  either  with  the  analysis  column 
or  in  taking  densities. 

It  is  now  desired  to  obtain  a  specific  value  of  PV 
at  60°C.  This  can  be  done  in  one  of  three  ways. 

(1)  The  variation  of  pressure  with  temperature  on  a 
constant  volume  of  gas  may  be  used  in  conjunction  with 
the  accepted  value  of  (PV)2^*  This  method  is  not 
considered  accurate  enough  for  the  present  measurements  a) 

(2)  The  value  of  (PoV0)^q  may  be  determined  directly  from 
(P0V0)2^  if  the  latter  is  known.  Since  this  procedure 

is  used  in  reverse,  to  determine  (P0V0)2j5  from  (PoVo)60 
in  this  work,  this  method  is  excluded. 

(3)  The  slope  of  the  line  (PV)/(PiV^)2^  vs  P  can  be  used 
in  conjunction  with  the  accepted  value  of  (P^Y-^)^  By 
this  method,  a  value  of  (Po^o)^q  -  1373*84  cc.atm./g. 

was  obtained.  Considering  the  errors  in  density  measurement, 
this  value  is  probably  accurate  to  closer  than  1  cc. 

The  experimental  values  of  (PV) ^o/(Pl^l)23  vs  P  are 
shown  in  Figure  2.  From  oaleulated  values  of  (P0V0)2^ 

and  (P o^Q)rjrj  and,  in  the  case  of  the  23°curve,  from  the 
calculated  slope,  tentative  compressibility  curves  at 
23a  and  77°  are  also  presented.  Data  obtained  by  Boyer 
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Figure  2 

Pressure  ~  Volume  -  Temperature 
Relationships. 
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and  Morrison  (1)  and  Gillies  and  Hugill  (6)  are  presented 
for  comparison. 

The  general  trend  of  Pv  vs  P  curves  can  he  predicted 
from  an  equation  of  state  such  as  that  of  Van  der  Waals : 

(P  *  A/V2)(V  -  B)  z  RT. 

Rearranged,  this  becomes  - 

PV  :  RT  «■  PB  -  A/V  «-  AB/V2 

At  low  pressures,  P  is  small  and  V  is  relatively  large. 
Hence  the  terra  a/v  is  considerably  more  significant  than 
(PB  ♦*  AB/V2).  Thus,  as  P  increases,  the  negative  value  of 
A/V  causes  PV  to  decrease.  However,  simultaneously  PB  is 
increasing  and,  due  to  the  decrease  in  V,  AB/V2  is  also 
increasing.  Due  to  the  fact  that  AB/V2  increases  at  a 
faster  rate  than  does  A/V,  the  positive  value  of  PB  ►  AB/V2 
soon  predominates  over  the  negative  value  of  A/V  and  the 
slope  of  the  curve  becomes  positive.  For  this  gas,  this 
occurs  in  the  regioxi  of  1.50  atmospheres. 

By  comparison  of  the  compressibility  curves  for  the 
three  gases  at  60°C,  it  will  be  evident  that  the  eas^of 
compressibility  increases  from  Royalite  #27  to  Alberta 
Oil  Incomes  #2  to  N. W.  Hudson’s  Bay  #6  in  the  low 
pressure  region.  This  may  be  explained  in  view  of  the 
different  compositions  of  the  three  gases. 

It  will  be  noted  that  an  increase  in  the  content  of 
propane  and  pentane  occurs  from  Royalite  #27  to  Alberta 
Oil  Incomes  #2  to  N.W.  Hudson^  Bay  #6.  Royalite  #27  and 
Alberta  Oil  Incomes  #2  have  the  same  butane  content  but 
that  of  U.  <V.  Hudson’s  Bay  #6  is  greater.  It  is  thus 
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evident  that  the  oonoentration  of  higher  hydrooarbons 
inoreases  progressively  from  Royalite  ft 27  to  Alberta  Oil 
Incomes  ftZ  to  N. W.  Hudson*s  Bay  ft6.  The  average 
concentration  as  shown  by  the  molecular  weight  of  the  gas 
shows  the  same  regular  trend. 


Molecular  weight 

Density 

(from  analyses) 

g/litre 

Royalite  ftZl • 

19.  96 

0.8006 

Alberta  Oil  Incomes  ftZ 

20.33 

0.8278 

N. y. Hudson*  s  Bay  ft6. 

20.49 

0.8131 

The  gas  densities,  shown  for  comparison,  are  not 
<iuite  so  regular,  but  the  low  temperature  fractionat ion 
is  considered  to  give  a  more  ideal  result. 


The  trend  of  the  PV  vs  P  curve  is,  as  has  been 
discussed  above,  the  result  of  two  opposing  factors: 

1.  The  force  of  intermolecular  attraction  tending  to 
decrease  the  value  of  PV  due  to  the  negative  sign  of  the 
term  A/V. 

2.  The  volume  occupied  by  the  molecules  tending  to 
increase  the  value  of  PV  due  to  the  positive  sign  of  the 
term  PB.  The  relative  values  of  A  and  B  will  determine 
the  effect  of  the  term  AB/V2 • 

In  the  case  of  the  gases  under  consideration,  it 
might  be  expected  that  the  increase  in  the  content  of  higher 
hydrocarbons  from  Royalite  to  U. W.  Hudson* s  Bay  ft& 
would  tend  to  increase  the  effect  of  the  ~A/V  term.  In 
regions  of  moderately  low  pressure,  where  the  volume  of 
the  molecules  is  relatively  unimportant,  the  inter¬ 
molecular  attraction  between  large  molecules  such  as 
pentane  would  be  expected  to  be  much  greater  than  similar 
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forces  between  the  smaller  methane  molecules.  Hence, 

increasing  the  pentane  content  of  the  gas  would  increase 

the  importance  of  the  term  involving  intermolecular 

attraction.  Provided  the  pressure  was  moderately  low, 

the  effect  of  the  larger  pentane  molecule  on  the  terra 

involving  molecular  volume  would  be  insignificant.  Thus, 

at  low  pressures,  the  reason  for  the  relative  positions 

of  the  three  PV  vs  P  curves  is  evident. 

As  pressure  increases,  the  molecular  volume  term 

becomes  relatively  more  important,  and,  as  has  been  seen 

before,  a  pressure  is  reached  at  which  (PB  *  AB/v2) 

increases  more  rapidly  than  does  A/V  and  a  minimum  in  the 

PV  curve  is  reached.  In  a  gas  having  a  higher  pentane 

content,  this  pressure  would  be  reached  sooner,  due  to 

the  effect  of  the  larger  pentane  molecules.  This  is 

evident  in  the  curves  illustrated.  The  line  for  H. W. 

Hudson1 s  Bay  #6  reaches  a  minimum  and  has  started  to 

increase  before  the  minimum  is  reached  in  either  the 

Alberta  Oil  Incomes  jf 2  curve  or  the  Royalite  #27  curve. 

The  effect  of  decreasing  temperature  is  to  increase 

the  relative  importance  of  the  molecular  attraction  term. 

This  is  evidenced  by  the  increase  in  the  negative  slope 

of  the  PV  curve  in  going  from  77°  to  60°  to  23°*  The 

o 

general  trend  of  the  three  curves  at  23  is  the  same. 

The  position  of  the  curve  at  77°  for  H.W.  Hudson* s  Bay  #6 

is  in  agreement  with  what  might  be  expected  from  the 

o 

behavior  of  Alberta  Oil  Incomes  #2  at  98  as  measured 
by  Boyer  and  Morrison(l) 
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Reference  was  made  earlier  to  the  deviation  per  om. 
obtained  by  Boyer  and  Morrison.  This  was  obtained  from 
the  expression  - 

P0V0/PiV1  =  1  «•  A  :  1  *•  bflj 
Where 

A  =  deviation  factor 
b  z  deviation  factor  per  cm. 

The  value  of  PV  at  any  pressure  and  temperature  within 
the  experimental  range  is  readily  obtained  from  large 
scale  plots  of  (PV)^/(P^ V^)2^  vs  P.  Knowing  the  slope  of 
the  curve  at  the  desired  point,  it  is  only  necessary  to 
multiply  by  the  accepted  value  of  (P^V^)^  to  get  an 
absolute  value  of  (PV)T.  This,  in  turn,  yields  a  value 
of  V,  ie.  the  specific  volume,  which  may  be  inverted  to 
give  the  density. 

In  addition  to  this  direct  use  of  the  PV  data,  they 
are  of  value  in  fomulating  equations  of  state,  the 
coefficients  of  which  lead  to  interatomic  forces  and 
theoretical  considerations  of  the  gaseous  state.  It  is 
also  possible,  although  rather  difficult,  to  calculate 
thermodynamic  properties  of  gases  from  the  pressure- 
volume-temperature  data.  A  detailed  description  of  the 
method  involved  is  given  by  Hewitt  (14).  This,  in  most 
cases,  requires  a  knowledge  of  V  -  T  coefficients. 

Graphs  of  V  vs  T  yield  very  inaccurate  values  of  the 
coefficients  and  it  is  highly  preferable  to  use  the 
variation  in  ‘'residual  volume"  rather  than  actual  volume. 
This  residual  volume  has  already  been  discussed  in  more 
detail  previously.  Using  the  method  outlined  above,  the 
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residual  volumes  were  calculated  for  N. W.  Hudson’s  Bay  #6 
✓  o 

gas  at  60  C.  The  results  are  shown  in  Figure  3  and 
Table  4.  For  comparison,  the  figures  for  Alberta  Oil 
Incomes  #2  and  Royalite  #27  are  included. 

It  will  be  noticed  that  in  all  three  cases,  "a"  is 
positive  at  low  pressures  and  follows  a  general  trend  to 
negative  values  at  high  pressures.  The  residual  volume, 

"a",  will  necessarily  be  zero  when  P0V0  =  PV. 

The  trend  of  "a",  from  positive  to  negative  values 
can  likewise  be  predicted  from  consideration  of  a  suitable 
equation  of  state  such  as  that  of  Van  der  Waals : 

(P  «■  A/V2)(V  -  B)  r  ET. 

Rearranging,  we  get  a  =  RT/P  -  V  =  A /PV  -  AB/PV2  -  B. 
When  pressure  is  low,  V  is  large  and  the  term  AB/PV2  can 
be  neglected.  Also  PV  changes  only  slightly  with  small 
changes  in  pressure.  Hence 

a  z  A/PV  -  B  =  Constant. 

Actually  PV  does  change  for  small  changes  in  pressure, 
but  small  variations  in  A  and  B,  compensate  and  keep  Mart 
constant.  This  is  readily  seen  in  the  graph  of 
a  vs  P  in  the  region  0-40  atmospheres.  As  pressure  is 
increased,  V  becomes  smaller  and  the  term  AB/PV2  soon 
predominates  over  A/PV  due  to  the  second  power  of  V. 

Thus  a  decrease  in  "a"  can  be  predicted  with  increase  in 
pressure  and  ultimately  ’’a11  will  become  negative. 

Changes  in  "a”  with  temperature  can  be  predicted 
roughly  by  substituting  RT  for  PV  in  the  Van  der  Waal’s 
expression.  Thus 


a  =  A/RT  "  ABP/R2T2  -  B 
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P.  a 

1 

20 

40 

60 

80 

100 

120 

140 

160 

180 

200 


Table  4 


Residual  Volumes,  oo/g,  at  60°C 


N.W.H.B./6 

2.623 

2.62? 

2.629 

2.378 

2.472 

2.296 

2.110 

1.923 

1.710 

1.498 


a.o.i.  n 


2*338 

2.311 

2.433 

2.341 

2.193 

2.048 

1.878 

1.703 

1.314 

I.303 


Royalite  27 

2.449 

2.449 

2.426 

2.344 

2.238 

2.113 

1.964 

1.780 

1.990 

1.397 

1.202 


•DUAL  VOLUME 


Figure  3 

The  Variation  of  the  " Residual  Volume" 


with  Pressure. 
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An  increase  in  temperature  will  obviously  decrease  "aM 
and  a  point  will  eventually  be  reaohed  where  "aM  is  no 
longer  positive  for  any  value  of  pressure.  The  point  at 
which  a  =  0  when  P  *  0  is  called  the  Boyle  Point. 

Changes  in  residual  volume  with  composition  can  be 
seen  by  a  similar  treatment  to  that  used  above  for  the 
compressibility  curves.  In  the  low  pressure  region, 
where  changes  in  B  are  negligible  for  increase  in  the 
proportion  of  higher  hydrocarbons  in  the  gas,  the 
increasing  importance  of  A  in  going  from  Royalite  -Jj2  7 
gas  to  H.W.  Hudson* s  Bay  ^6  gas  will  cause  the  residual 
volume  to  increase.  This  is  evidenced  by  the  higher 
initial  value  of  '‘a"  for  the  latter  gas.  As  the  pressure 
increases,  the  larger  concentration  of  higher  hydrocarbons 
in  H.W.  Hudson* s  Bay  ^6  gas  causes  the  effect  of 
molecular  volume  to  be  felt  at  a  lower  pressure  than  is 
the  case  for  the  other  two  gases.  Thus  it  is  not 
surprising  that  the  curve  for  H.W.  Hudson* s  Bay  ^6  gas 
crosses  the  curve  for  Alberta  Oil  Incomes  jj 2  gas  at  164 
atmospheres.  At  a  still  higher  pressure  it  is  to  be 
expected  that  txhis  curve  will  likewise  cross  the  curve  for 
Royal  it e  ^27 

The  preceding  discussions  on  the  behaviour  of  the  gas 
as  predicted  by  Van  der  Waals  is  not  quantitative 
because  of  the  limited  accuracy  of  the  Van  der  .Vaals* 
equation  of  state.  A  more  exact  analysis  would  necessitate 
a  much  more  complex  equation  of  state  and  the  difficulties 
involved  nave  already  been  pointed  out.  It  has  also  been 
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shown  that  a  quantity  called  fugacity  can  be  employed  to 
render  the  simple  equations  exact  for  a  real  gas. 

Reference  has  also  been  made  to  an  easy  method  by  which 
fugacities  c-.n  be  calculated  from  residual  volumes 
determined  from  pressure-volurae-temperature  measurements. 

Fugacities  calculated  for  N. W.  Hudson* s  Bay  #6  gas 
are  shown  in  Table  and  illustrated  graphically  in 
Figure  4.  Fugacities  for  Alberta  Oil  Incomes  ft 2  are 
shown  for  comparison,  in  Table  5. 

From  the  effect  of  composition  on  the  residual  volume, 
it  would  be  expected  that  the  larger  proportion  of  higher 
hydrocarbons  in  the  N.W.  Hudson* s  Bay  ftb  gas  would  result 
in  a  lower  measured  fugacity  at  the  lower  pressures. 

This  is  readily  seen  to  be  the  case.  As  the  pressure 
increases,  the  effect  of  the  molecular  volume  term,  B, 
in  the  Van  der  Waals  expression  becomes  apparent  and  the 
residual  volume  for  IT. W.  Hudson* s  Bay  ft(>  decreases  much 
more  rapidly  than  that  for  Alberta  Oil  Incomes  ftZ*  A 
pressure  will  eventually  be  reached  at  which  the  areas 
under  the  two  a  vs  P  curves  are  the  same  and  the 
measured  fugacities  at  this  pressure  will  be  the  same  for 
the  two  gases.  With  further  increase  in  pressure,  the 
fugacity  of  IT.  w.  Hudson*  s  Bay  ftG  gas  will  exceed  that  of 
Alberta  Oil  Incomes  ftZ  gas. 

The  absolute  differences  between  the  fugacities  of 
the  two  gases  at  any  pressure  within  the  range  studied, 
is  slight  as  would  be  expected  in  view  of  their  close 
similarity  in  composition. 
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Table  3 


Fugacities  at  60°C 


P.  atm.  N.W.  Hudson* s  Bay  $6.  Alberta  Oil  Incomes^ 


fl 

f2 

fl 

f2 

1 

0,998 

20 

19,23 

18.79 

19*26 

18.86 

40 

37.03 

35.70 

37.10 

33.87 

60 

33.31 

51.55 

33.62 

31.74 

80 

68.77 

66.40 

68.98 

66 « 60 

100 

83.03 

8O.63 

83.36 

80.79 

120 

96.48 

94.43 

96.92 

94.36 

140 

109*30 

107.91 

109.82 

107.96 

160 

121.63 

121.14 

122.21 

121.14 

180 

133*70 

134.22 

134.23 

134.18 

200 

146.04 

147.18 

f^  =  experimental  fugacity 

values 

f2  z  fugacities 

calculated 

from  fugacities 

of 

pure  components  using 

Lewis*  rule. 

FUGACITV  -  ATMOSPHEK 


Figure  4 

The  Variation  of  Fugacity  with 


Pressure. 
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The  fugaoities  of  the  pure  components,  used  in  the 
application  of  Lewis1  rule,  are  the  same  as  those  used 
by  Boyer  and  Morrison  and  are,  in  the  main,  from  the 
investigations  of  Sage  and  Lacey.  As  shown  by  Boyer,  the 
application  of  Lewis1  rule  in  calculation  of  fugacities 
of  gaseous  mixtures  is  of  considerable  value,  especially 
in  engineering  calculations. 


Symbols 


Wg  -  Grams  of  gas  as  butanes  and  lighter. 

WD  "  Total  grams  in  phase. 

(Vg)-  Volume  of  gas,  butanes  and  lighter,  at  t°C  and  1  atm. 

(VQ)-  Volume  of  oil,  pentanes  and  heavier,  in  same  units 
as  (Vg). 

Fg  *  Cubic  feet  of  gas,  butanes  and  ligher,  at  20°C,  1  atm. 

Bbl-  33  Imperial  gallons. 

WQ  -  G.  of  oil,  including  pentanes  and  heavier. 

Vp  -  Volume  of  phase  under  equilibrium  conditions, 

(Vp)g  for  gas  phase;  (Vp)j,  for  liquid  phase. 

GQ  -  Gallons  of  oil,  pentanes  and  heavier,  at  20°C. 

Fp  -  Volume  of  phase  under  equilibrium  conditions  in 
cubic  feet. 

o 

Xp  -  Volume  of  phase  when  expanded  to  1  atm.  and  23  C. 

-  Ratio  of  grams  of  methane  per  cc.  of  gas  phase  to 
grams  of  methane  per  cc.  of  liquid  phase. 

Kg  ,Kz  jKjj  similar  ratios  for  ethane,  propane,  butane 
pand  nitrogen. 

Kc  ~  Ratio  of  grams  of  light  fraction  (CcH^p*)  per  cc. 
of  gas  phase  to  grams  of  same  fraction^as 
determined  by  refractive  index)  per  cc.  of  liquid 
phase. 

Kl"n  Ratio  of  grams  of  methane  nitrogen  per  cc.  of  gas 
phase  to  grams  of  methane  *  nitrogen  per  eg.  of 
liquid  phase. 


Table  6 


General  Experimental  Data 

70°  77°  85° 

Liquid  Phase  Samples:  59*59  co.  59*61  co.  59*62  cc. 
Gas  Phase  samples:  27*94  cc.  27*95  co*  27*97  cc. 


Original  Materials  (as  received  in  the  laboratory) 


Run  No .  1 

Pressure  at  equilibrium,  atm.  228.1 

Temperature  at  equilibrium  76.4 

Oil,  cc.  at  25°C  692.55 

Oil,  g.  •  561.71 

Gas,  cc.  at  25°C,  1  atm.  151996 

Gas,  g-  Q  .  0  125.89 


cc  gas,  25°C,  1  atm/cc  oil,  25  C0  219*5 
Cu.  ft.  gas,25  C.latm/bbl  oil  25  C  1218. 9 


2 

1F7T5 

77*0 

596.52 

483.84 

137915 


112.41 

231.2 

1283.9 


152T8 

77*0 

621.25 

503.90 

145559 

118.47 

255*9 

1298.9 


5 

2cyn 

77.0 

656.52 

552.51 
155998 

125.52 


234.6 

1302.7 


Original  Materials  Charged. 


Oil,  cc 

0 

pentane  ►,  25 

697.16 

Oil,  g. 

564.59 

Gas,  cc 

Q 

butane  -,  25  ,  1  atm. 

151076 

Gas,  g. 

121.01 

cc.  butane  -  25  »  1  atm/  oc 

216.7 

pentane  *,  25 
butane  -  /bbl  pentane  * 

Cu.  ft. 

I203.4 

600.70 

486.45 

137081 

109.80 

228.2 

1267*2 


625.69  661.22 
506.66  559*43 
144460  153066 
115.71  122.60 
20C.  9  231.4 

1282.2  1285.0 


6 

iwri 

70.0 

614.44 

498.38 

141164 

115.06 

229*7 

1275.6 


618.75 

501.06 

140310 

112.38 

226.8 

1259.5 


7 

1BTJT0 

85.0 

615.10 

498.91 

139560 

113.76 

226.9 

1260.0 


619.38 

501.57 

138716 

111.10 

224.0 


1243.9 


Liquid  Phase  Sample 
Liquid  Separated 


«•  ,  0 
density,  g/co,  20  C 

volume,  cc,  20°C 
molecular  weight 
volume,  light  fraction,  co 
weight,  light  fraction,  g. 

Gas  volume,  oc,  25°C,  1  a$m. 

Gas  density,  g/litre,  25  0,  1  atm 
Gas  mass,  g. 

Composition  of  gas,  mole  t 

30.58 

0.8298 

56.85 

178.0 

4.75 

6831.51 

0.8625 

5.8922 

£H4 

75.77 

Ng 

2.91 

c?h6 

12.22 

CzHg 

6.25 

C4%0 

C?Bf2». 

1.89 

0.96 

Gas  Phase  Sample 

Liquid  Separated 

g*  o 

0.5728 

refractive  index,  20  C 

1.3969 

molecular  weight 

94.5 

density,  g/cc.  20°C 

0.7107 

Gas 

volume,  cc,  25°C,  1  atm. 

6329.99 

Gas 

density,  g/litre,  25°C,  1  atm. 

0.7876 

Gas 

mass,  g. 

4.9855 

Composition  of  gas,  mole  % 

ch4 

c|h6 

& 

C5H12  ► 


84.71 

5.95 

7.13 

2.50 

1.48 

0.65 


Liquid  Phase  Sample 


Liquid  Separated 

1 

2 

g.  o 

30.58 

35.90 

density,  g/cc,  20  C 

0.82?8 

0.8188 

volume,  cc,  20°C 

36.85 

43.84 

molecular  weight 

178.0 

187.5 

volume,  light  fraction,  cc 

4.75 

6.09 

weight,  light  fraction,  g. 

3.38 

4.14 

Gas  volume,  cc,  25°C,  1  a£m. 

6831.51 

3320.75 

Gas  density,  g/litre,  25  0,  1  atm 

0.8625 

0.8899 

Gas  mass,  g. 

5.8922 

4  .  75  49 

Composition  of  gas,  mole  f. 

ch4 

75.77 

74.00 

n2 

2.91 

1.66 

CpH/- 

12.22 

14.07 

c,h§ 

6.25 

6.59 

C4%0 

C^HTp*- 

1.89 

0.96 

2.26 

1.42 

Gas  Phase  Sample 

Liquid  Separated 

g*  0 

0.5728 

O.3324 

refractive  index,  20  C 

1.3969 

1.3931 

molecular  weight 

94.5 

92.1 

density,  g/cc.  20°C 

0.7107 

0.6798 

Gas  volume,  cc,  25°C,  1  atm. 

6329.99 

4713.09 

Gas  density,  g/litre,  25°C,  1  atm. 

0.7876 

0.7669 

Gas  mass,  g. 

4.9855 

3.6145 

Composition  of  gas,  mole  % 

ch4 

84.71 

85.16 

u2 

3.53 

3.64 

c|h6 

7.13 

6.90 

C4HIO 

2.50 

2.40 

1.48 

1.11 

C5H12  ► 

O.65 

0.79 

-J- 

6 

7 

34.68 

34.20 

36.92 

36.50 

0.8225 

0.8205 

0.8211 

0.8227 

42.16 

41.68 

44.96 

44.37 

186.1 

188.0 

175.0 

is* 

183.2 

6.37 

5.08 

7.42 

4.48 

4.49 

5879.68 

3.48 

6248.48 

5.21 

5552.17 

3.23 

5462.12 

0.9477 

0.8884 

0.9300 

0.9216 

5.5722 

5.5511 

5.1635 

5.0339 

71.11 

73*46 

71.11 

72.08 

2. 66 

3.30 

3*12 

3.49 

13-23 

12.42 

13.44 

13.01 

6.84 

6.90 

7.34 

6.87 

4.23 

2.58 

3.45 

3.08 

1.93 

1.34 

1.54 

1.47 

0.4956 

0.5670 

0.3851 

0.3202 

1.3982 

1.3936 

1.3961 

1.3869 

95.4 

92.4 

94.0 

88.8 

O.705O 

0. 6856 

O.7015 

0.7212 

5066.75 

5611.01 

4826.97 

4777.30 

0.7442 

0.7608 

0.7399 

0.7636 

3.7707 

4.2689 

3.5715 

3.6*79 

84.17 

85.36 

85.86 

85.10 

4.99 

4.04 

4.18 

4.01 

7.18 

7.04 

6.82 

7.14 

2.34 

2.29 

2.14 

2.22 

0.93 

0.77 

0.  64 

0.93 

0.39 

0.50 

O.36 

0.60 

Table  7 


Properties  of  Gas  and  Liquid  Phases  and 
of  the  System  at  Equilibrium. 


Run  No. 

1 

2 

3 

5 

6 

Pressure,  atm. 

278*.  1 

rsr.5 

18278 

?5J*l 

157. 7 

Temperature. 

76.4 

77.0 

77.0 

77.0 

70.0 

Liquid  Phase 

Volume,  cc  (V_.)T 

1070.8 

783.4 

833.5 

' 

903.6 

779.5 

Mass,  g.  P  L 

655.12 

534.OO 

562.78 

602.61 

550.48 

Density,  g/cc. 

0.6118 

0.6817 

0.6752 

u . 6669 

0.70&2 

Composition  in  Weight  % 

1 

CH 

9.26 

6.34 

6.86 

7.49 

6.18 

N2 

0.62 

O.25 

O.45 

0.59 

0.47 

Cpll6 

2.80 

2.26 

2.40 

2.38 

2.19 

C^Eq 

2.10 

1.55 

1.82 

' 

1.94 

1.75 

°4H10 

0.84 

0.70 

1.48 

. 

• 

0.95 

1.08 

Lignt  fraction 

9.80 

10.74 

11.99 

9.37 

12.98 

Residue 

74.58 

78.16 

75.00 

77.28 

75.35 

w_/w_ 

O.I563 

0.1110 

0.I300 

0. 3.335 

0.1167 

)fn/( V  )T 

111.60 

86.52 

95.11 

101.68 

90.20 

(vg)f§/(vpn;0 

179.07 

116.68 

132.78 

144.08 

118.50 

994.40 

647.94 

737.35 

j 

800.10 

658.05 

Gas  Phase 

Volume  cc,  (V  ) 

153.2, 

440. 6 

390.5 

320.4 

444.5 

Mass,  g.  p  s 

30.47 

62.22 

59.59 

55.43 

62.94 

Density,  g/cc. 

0.1989 

0.1412 

0.1526 

0.1730 

0.1416 

Composition  in  Weight  % 


fraction 


in  cc. 


G0/iL  iS  1000  cu.  ft. 
ip/(Yp)g  in  same  units. 


63.63 
4.  63 
10.04 
5*17 
4.O3 
12.30 

0.0249 

0.0037 

0.0244 

231.82 


63.73 

4.91 

9.?8 

5.09 

3.10 

11.17 

0.0138 

0.0038 

0.0220 

171.78 


63.IO 

6.33 
10.08 

4.83 

2.33 
12.93 


0.0197 

0.0043 

0.0243 

I83.83 


64.44 

5*32 

9.96 

4.75 

2.11 

13.42 

0.0232 

0.0034 

0.0264 

206.09 


66.9] 
5*  7< 
9.9' 

4.3! 
1.8< 
11. 0( 

O.OI31 

O.OO31 

0.020- 

176.3] 


7 

I8T7.0 

83.0 


791.4 

531.29 

0.6966 


6.28 


0.53 

2.12 

1.64 

0.97 


8.33 

80.11 


O.H33 

88.76 

II8.23 

656.66 


System  Properties 


(Vp)g/(Vp)1 

Equilibrium 

l 

07IT31 

Constants 

2 

O.5625 

07^685 

5 

0  *  354 6 

6 

075702 

7 

0. 5466 

1 

2.233 

2.148 

2.075 

2.230 

2.174 

2.158 

1.170 

0.916 

0.9J7 

I.089 

0.916 

1.021 

0.798 

0.679 

0.602 

0.636 

0.524 

0.596 

M 

1.569 

0.917 

O.390 

0.578 

0.338 

0.544 

£c 

0.416 

0.216 

0.243 

0.371 

0.170 

0.249 

% 

2.421 

4.059 

3.333 

2.359 

2.455 

2.108 

Kl+Jf 

2.244 

2.220 

2.151 

2.239 

2.194 

2.154 

432.6 

61.39 

0.1419 


66.49 

5.47 

10.46 

4.77 

2.62 

10.18 

0.0145 

0.0033 

0.0190 

173.95 


903. 6 
602.61 
0.6669 


7.49 

0.^9 


2.38 

1.94 

0.95 

9-37 

77.28 


0 •1335 
101.68 
144.08 
800.10 


320.4 

55*43 

0.1730 


64.44 

5.52 

9.96 

4.75 

2.11 

13.42 


0.0232 

0.0054 

0.0264 

206.09 
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Equilibrium  Measurements 

The  direct  experimentally-determined  data  are 
presented  in  Table  b.  The  derived  data  are  shown  in 
Table  7*  In  experiments  1,  2,3  and  3  an  attempt  was 
made  to  approach  the  reservoir  conditions  of  temperature 
and  gas: oil  ratio,  while  allowing  the  equilibrium  pressure 
to  vary  above  and  below  that  existing  in  the  reservoir. 

In  experiments  6  and  7  an  attempt  was  made  to  maintain 
the  reservoir  conditions  of  pressure  and  gas  roil  ratio 
while  allowing  temperature  to  vary  above  and  below  that 
existing  in  the  reservoir.  A  further  series  of  two 
experiments  was  done  in  which  an  attempt  was  made  to 
maintain  the  reservoir  conditions  of  pressure  and  temp¬ 
erature  while  allowing  the  gas  roil  ratio  to  vary. 
Unfortunately  no  results  were  obtained  in  this  last  series 
due  to  failure  of  the  high  pressure  apparatus  during  the 
analyses.  Temperature  control  was  very  exact  in  all 
cases  but  it  was  much  more  difficult  to  introduce  the 
size  of  charge  which  would  give  simultaneously  a 
desired  pressure  and  a  desired  gas  roil  ratio.  In  all 
cases  except  in  the  first  table  on  '’original  materials 
(as  received  in  the  laboratory)”  the  term  ”gas”  refers 
to  butanes  and  lighter,  and  ”oil”  refers  to  pentanes  and 
heavier  const ituents. 

The  weight  of  liquid  separated  from  the  liquid  phase 
was  found  to  fc0.01  g.  and  its  density  was  determined 
gravimetric ally .  The  molecular  weight  of  the  oil  was 
determined  by  the  cryoscopic  method  but  so  much  variation 
was  found  in  the  results  that  it  was  considered 
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Inadvisable  to  use  the  raoleoular  weights  so  determined. 

The  reason  for  this  variation  is  still  obscure  but  may 
be  due  to  traoes  of  water  condensing  in  parts  of  the 
freezing  point  apparatus  and  thereby  changing  the 
temperature  gradient  by  impeding  the  rate  of  heat  transfer. 

The  molecular  weight  of  the  light  fraotion  from  the 
gas  phase  sample  is  considered  to  be  fairly  accurate. 

The  reason  for  this  is  that  the  molecular  weight  was 
determined  for  each  sample  by  comparing  its  refractive 
index  with  a  refractive  index  vs  molecular  weight  graph 
previously  prepared.  This  graph  was  prepared  by  measuring 
molecular  weights  of  different  fractions  from  a  series 
of  A.S.T.M.  distillations  and  represents  the  best 
straight  line  through  some  lj?  points  in  a  fairly  narrow 
molecular  weight  range  (Figure  1).  However  in  view  of 
the  unreliability  of  the  molecular  weights  of  liquid 
from  the  liquid  phase  it  was  considered  advisable  to  deal 
only  with  weight  fractions  for  each  phase,  rather  than 
mole  fractions. 

The  weight  of  the  light  fraction  from  the  gas  phase 
was  accurately  determined  by  weighing  the  sealed 
dilatometer  and  its  volume  was  determined  by  immersing 
the  dilatometer  in  water  at  20°C  and  reading  the  volume 
directly  from  calibrations  etched  on  the  dilatometer. 

The  corresponding  light  fraction  in  the  liquid  from  the 
liquid  phase  was  determined  by  fractional  distillation 
of  the  liquid  and  plotting  %  over  vs  refractive  index. 

The  volume  having  the  same  refractive  index  as  the 
dilatometer  sample  was  considered  to  be  light  fraction. 
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All  gas  volumes  were  measured  in  the  vicinity  of 

jj.3-78  om.  of  Moroury  at  room  temperature  and  oorrected 

using’  the  deviation  faotor  for  N. «/.  Hudson»s  Bay  #6  gas. 

Since  it  is  not  possible  to  determine  mole  fractions 

on  the  phase  samples,  the  usual  equilibrium  constants 

are  replaced  by  constants  determined  from  concentrations. 

Thus,  in  this  present  case, 

Kf  z  grams  methane  per  cc.  of  gas  phase 

grams  methane  per  cc.  of  liquid  phase. 

This  K  is  analogous  to  the  constant  in  Henry1 s  Law  and 

becomes  Kjg  on  multiplication  by  the  ratio  of  molar  volumes 

in  the  corresponding  phases. 

Considering  the  data  in  Table  7,  we  see  that  the 

solubility  of  butanes  and  lighter  in  pentanes  and 

heavier  components,  expressed  as  cubic  feet  per  barrel 

(Fg/bbl  increases  with  increasing  pressure  and 

the  rate  of  increase  is  greater  at  the  higher  pressures. 

This  is  what  would  be  expected  from  a  consideration  of 

solubility  work  done  by  other  workers. 

Considering  runs  ^2,  6  and  2,  7  it  is  evident  that 

increasing  temperature  decreases  the  solubility  of  butanes 

and  lighter  in  pentanes  and  heavier  constituents. 

770  -  8.5°  than  from 

However  the  decrease  is  larger  from^O0  “  77°.  This  may 
be  partially  explained  by  consideration  of  the  other 
variables  involved.  Three  factors  simultaneously  affect 
the  solubility  values.  Consider,  first,  runs  #2  and  ^6. 

(1)  There  is  a  decrease  in  pressure  of  about  2.2 
atmospheres  in  going  from  ,70°  -  77°  which  will  cause  a 
slight  decrease  in  solubility.  (2)  The  gas:oil  ratio 
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increases  slightly.  This  will  have  the  effect  of  increasing 
the  volume  oooupied  by  the  gas  phase,  and  henoe,  there 
will  be  a  more  complete  transfer  of  the  lighter  liquid 
hydrocarbons  to  the  gas  phase.  This  will  leave  a  heavier 
liquid  residue  and  the  butanes  and  lighter  components 
will  have  a  decreased  solubility.  (3)  The  increase  in 
temperature  will  have  two  effects :  the  normal  effect  of 
temperature  on  solubility  —  at  low  pressures,  a 
temperature  increase  will  cause  a  decrease  in  solubility; 
as  the  pressure  increases  a  region  will  be  found  in  which 
the  solubility  decreases  to  a  minimum  and  then  increases 
with  increasing  temperature.  The  other  effect  of 
increasing  temperature  is  due  to  the  increased  vapor 
pressure  of  the  light  liquid  fractions.  Vaporization  of 
these  components  leaves  a  heavier  liquid  residue  in  which 
the  gaseous  compbnents  have  a  decreased  solubility. 

Consider,  now,  runs  #3  and  ^7*  The  decrease  in 
pressure  in  going  from  77°  to  83°  is  2*8  atmospheres 
which  is,  for  practical  purposes,  identical  with  the 

pressure  effect  on  solubility  in  going  from  70°  to  77°. 

0 

The  gas:oil  ratio,  however,  decreases  in  going  from  77 
to  83°  so  the  effect  will  be  to  cause  an  increase  in  gas 
solubility.  However,  in  both  temperature  ranges  the 
effect  of  such  a  small  change  in  gas:oil  ratio  will  be 
very  slight.  In  the  higher  temperature  range  (77°to  83°) 
the  increase  in  vapor  pressure  of  light  liquid 
components  witn  increase  in  temperature  will  be 
somewhat  larger  than  in  the  lower  temperature  range. 

Thus,  one  might  expect  a  smaller  decrease  in  solubility 
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between  70°  •  77°  than  between  77°  “  83°.  A  further 
variable  will  be  the  absolute  value  of  the  pressure. 

As  observed  by  Johnson  (8),  at  low  pressures  the  solubility 
vs  temperature  ourve  has  only  a  slight  ourvature,  which 
becomes  much  more  pronounced  at  higher  pressures.  Thus, 
it  may  be  that  the  larger  decrease  in  solubility  between 
77°  "  83°  than  between  70°  -  77°,  is  partially  due  to 
the  fact  that  the  former  values  are  at  a  pressure  of 
180  atmospheres  as  compared  to  167  atmospheres  for  the 
latter  values. 

The  closest  comparable  data  for  Alberta  Oil  Incomes 
fZ  are  at  I38.6  atmospheres  and  98°C.  It  was  considered 
inadvisable  to  extrapolate  the  data  for  H.W.  Hudson’s  Bay 
j^6  over  such  a  wide  range  of  temperature  and  pressure. 

Considering  the  vaporization  of  pentanes  and  heavier 
constituents  into  the  gas  phase,  for  H.W.  Hudson’s  Bay  ^6, 
expressed  as  gallons  of  pentanes  and  heavier  per  thousand 
cubic  feet  of  expanded  phase  (G0/Xp),  it  is  apparent, 
from  Runs  2,  3  and  3,  that  an  increase  in  pressure  causes 
an  increase  in  the  pentane  *■  content  of  the  gas  phase. 

A  clear  distinction  should  be  made,  at  this  point, 
between  gallons  of  pentanes  and  heavier  per  cubic  foot 
of  phase  at  equilibrium  (G0/(Fp)g)  and  Per  cul3ic  f00*  of 


The  former  value  will  always 


show  an  increase  with  pressure  as  predicted  by 
Poynting.  This  is  evident  from  the  values  of  (G’0/$’p)g) 

in  Table  7»  However,  the  latter  value  (G0/Xp)  will  show 
an  increase  with  pressure  only  when  the  transfer  of 

light  liquid  components  to  the  gas  phase  occurs  at  a 
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faster  rate  than  oan  be  compensated  for  by  the  increased 

concentration  of  gaseous  components  in  the  gas  phase  due 

to  compression.  Thus  G0/Xp  will  increase  with  pressure 

only  in  the  region  of  retrograde  condensation.  At  all 

lower  pressures  Go/Xp  will  decrease.  The  fact  that 

Go/^p  increases  in  going  from  167«.5  atmospheres  to  20^.1 

atmospheres  indicates  that  the  region  of  retrograde 
a 

condensation  is  encountered  in  this  pressure  range. 

Run  j^l  introduces  an  anomalous  result,  in  showing  a 
decrease  in  G0/Xp  with  still  further  increase  in 
pressure.  Since  this  experiment  was  the  first  one  carried 
out,  the  only  apparent  explanation  is  that  the  value  of 
G0/Xp  is  in  error  due  to  inexperience  in  handling  the 
analytical  apparatus. 

Considering  the  temperature  effect,  it  is  evident 
that  three  factors  affect  the  transfer  of  pentane  and 
heavier  components  to  the  gas  phase.  In  the  temperature 
range  70°  -  77°  (Runs  ^6  and  2)  there  is  a  slight 
decrease  in  pressure  which  should  cause  a  slight  decrease 
in  G0/Xpj  the  simultaneous  increase  in  gas:oil  ratio 
should  also  decrease  the  pentane  content  of  the  gas  phase. 
Although  a  more  complete  transfer  of  light  liquid 
hydrocarbons  to  the  gas  phase  might  be  expected  to 
to  accompany  an  increase  in  gas: oil  ratio,  the  diluting 
effect  of  the  lighter  gaseous  components  would  result  in 
a  net  decrease  in  G0/Xp.  Both  the  pressure  effect  and 
the  gas: oil  ratio  effect  will  have  a  very  slight  influence 
on  the  system.  In  opposition  to  these  effects,  is  the 
effect  due  to  increase  in  temperature,  which  should 
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produce  a  decided  increase  in  G0/Xp  due  to  the  increased 
vapor  pressure  of  the  volatile  liquid  components.  The 
increase  in  pentane  and  heavier  constituents  in  the  gas 
phase  due  t o  an  increase  in  temperature  would  be 
expected  to  predominate  over  the  contrary  effect  of 
pressure  decrease  and  gas: oil  ratio  increase.  This  would 
account  for  the  higher  value  of  G0/Xp  at  77°  (Run  #2) 
than  at  70°  (Run  #6). 

Comparing  Runs  3 ,  7  in  which  the  temperature  rises 
from  77°  to  83°,  it  is  evident  that  the  decrease  in 
pressure  will  decrease  the  value  of  G0/Xp  by  an  amount 
very  similar  to  the  effect  predicted  in  the  Runs  at 
70°  and  77°.  The  effect  here,  will  be  equally  small. 

The  simultaneous  decrease  in  gas:oil  ratio  will  have  a 
slight  tendency  to  increase  the  transfer  of  pentanes  and 
heavier  constituents  to  the  gas  phase,  as  expressed  by 
G0/Xp*  However,  here  too,  the  effect  of  increased 
temperature  in  increasing  the  vapor  pressure  of  the 
volatile  liquid  components  would  be  expected  to  predominate 
over  all  other  factors,  and  an  increase  in  G0/Xp  might  be 
predicted.  The  fact  that  a  very  decided  decrease  is 
measured,  would  lead  to  the  conclusion  that  Run  ^7  is 
very  much  in  error. 

The  reason  for  this  discrepancy  is  not  apparent. 

Brrors  in  measuring  the  pentane  *■  fraction  of  the  gas  from 

the  gas  phase  may  be  as  high  as  20fo  but  this  is  only  a 

L 

small  portion  (about  13-2Qf0?\volume)  of  the  "light  fraction" 
most  of  which  is  retained  and  accurately  measured  in  the 
dilatometer.  The  error  appears  to  be  in  Run  ^7  rather 
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than  in  any  of  the  other  runs.  A  possible  explanation 
for  the  anomalies  observed  in  Run  #7  will  be  suggested 
later. 

Reference  might  be  made,  here,  to  the  difficulty 
encountered  by  Boyer  and  Morrison(l)  in  the  appearance 
in  the  gas  phase  sample  of  a  liquid  whose  presence  could 
not  be  accounted  for.  Further  investigations  showed 
that  this  liquid  was  carried  as  spray  or  as  liquid  through 
the  lines  and  trapped  in  the  gas  phase  pipette  due  to  its 
horizontal  position.  In  this  work  an  improved  design  of 
the  apparatus  has  been  used,  in  which  the  gas  phase 
pipette  occupies  a  vertical  position  in  the  oil  bath. 

This  has  eliminated  all  liquid  phase,  as  such,  from  the 
gas  phase  sample.  The  liquid  which  does  separate  out  in 
the  dry  ice-acetone  condensers  was  present  in  the  pipette 
as  a  gas. 

Gas  phase  densities  show  a  decided  increase  with 
pressure  as  shown  by  comparison  of  Runs  1,  2,3  s 

This  is  a  result  of  three  factors:  (1)  normal  increase 
in  density  due  to  compression  of  the  gas  with  increased 
pressure  (2)  preferential  solubility  of  propane  and 
butane  in  the  liquid  phase  from  the  gas  phase  (3)  transfer 
of  pentanes  and  heavier  components  from  the  liquid  phase 
to  the  gas  phase.  The  second  of  these  factors  tends  to 
deorease  the  gas  phase  density  in  opposition  to  the 
combined  effects  of  the  other  two  factors  tending  to 
increase  the  gas  phase  density.  The  first  and  third 
factors  predominate.  The  net  result  is  an  increase  in 
gas  phase  density  which  is  roughly  proportional  to  the 
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Variation  of  Phase  Density  with 


Pressure 
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increase  in  pressure.  The  more  rapid  rate  of  increase 
in  gas  phase  density  at  the  higher  pressure  may  be  due  to 
the  increased  rate  of  transfer  of  pentane  and  heavier 
constituents  from  the  liquid  phase. 

A  simultaneous  increase  in  gas  roil  ratio  might  be 
expected  to  decrease  the  gas  phase  density  due  to  the 
diluting  effect  of  the  gaseous  hydrocarbons  on  the 
volatile  liquid  components.  It  is  difficult  to  confirm 
this  suspected  trend  from  the  data,  due  to  the  small 
variation  in  gasroil  ratio  encountered.  The  effect  is 
too  slight  in  these  experiments  to  have  any  significance. 

The  temperature  effect  on  gas  phase  density  is 

shown  in  Runs  2,  6  and  3,  7*  There  is  a  two-fold  effect 

of  temperature  on  gas  phase  density:  (1)  at  constant 

pressure,  an  increase  in  temperature  would  cause  a  decrease 

in  density  which  is  roughly  proportional  to  the  change  in 

the  absolute  temperature.  This  is  a  direct  conclusion  from 

Charles*  law.  (2)  due  to  the  increased  vapor  pressure 

of  light  liquid  hydrocarbons,  an  increased  temperature 

should  increase  the  gas  phase  density.  If  the  phase 

density  at  77°  and  167*5  atmospheres  (Run  2)  is 

o 

calculated  from  that  at  70  and  169*7  atmospheres  (Run  6) 
assuming  a  linear  relationship  between  density  and 
pressure  and  correcting  for  temperature  on  the  basis  of 
the  change  in  the  absolute  temperature,  a  density  of 
0.1370  g/cc  is  calculated  for  Run  #2.  The  measured  value 
is  0.1412  g/cc.  This  difference  must  be  a  result  of  the 
increased  vapor  pressure  of  the  volatile  liquid 
hydrocarbons  at  the  higher  temperature,  tending  to  increase 
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the  gas  phase  density. 

Similar  calculation  of  the  gas  phase  density  in 

Run  #7  from  that  in  Run  yields  a  value  of  0.146?  g/oo. 
The  measured  value  is  0.141?  g/cc. 

It  is  reasonable  to  expect  tha't  the  calculated  phase 
density  at  83°C  would  be  lower  than  the  measured  value 
since  the  calculation  does  not  take  into  account  the 
increased  transfer  of  light  liquid  fractions  to  the 
gas  phase  at  the  higher  temperature.  The  reverse  is 
actually  the  case  here.  In  view  of  the  anomalies 
apparent  in  Run  §']  it  seems  reasonable  to  assume  that  the 
data  for  Run  fl  are  in  error. 

The  simple  assumption  that  a  state  of  equilibrium 
had  not  been  attained  at  the  time  of  taking  samples, 
serves  to  explain  all  the  discrepancies  which  have  been 
observed.  The  decrease  in  solubility  of  gas  in  oil,  as 
expressed  by  Fg/bbl  for  the  temperature  interval 

70°  -  77°  is  roughly  10  cu.  ft./bbl  as  compared  with  80 
cu.  ft./bbl  decrease  for  the  temperature  interval  77°"8^°. 
This  difference  is  easily  explained  if  the  oil  in  Run  f] 
had  not  become  saturated  with  gas.  Similarly,  in  the 
transfer  of  volatile  liquid  fractions  to  the  gas  phase 
with  increasing  temperature,  a  decrease  in  Go/%.  is 

.cr 

observed  between  77°  “  8^°  whereas  an  increase  had  been 
predicted  on  the  basis  of  the  behavior  of  the  systems  at 
70°  and  77°.  If  equilibrium  had  not  been  attained  in 
Run  #7,  the  transfer  of  volatile  liquid  fractions  to  the 
gas  phase  would  not  be  complete,  and  the  value  of  Go/^p 
would  be  lower  than  would  be  expected. 
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This  incomplete  transfer  of  pentane  and  heavier  components 
to  the  gas  phase  at  8^°C  due  to  non-attainment  of 
equilibrium,  would  also  explain  why  the  gas  phase  density 
at  85 °C  calculated  from  that  at  77°C,  is  higher  than  the 
experimental  density.  If  equilibrium  had  been  attained, 
it  is  reasonable  to  assume  that  the  measured  density  would 
have  been  much  higher  than  was  actually  the  case. 

The  anomalies  in  liquid  phase  density  to  be 
discussed  presently,  are  also  readily  accounted  for  by 
the  assumption  that  equilibrium  had  not  been  attained  in 
Run  j^7« 

It  is  possible  to  compare  Run  fZ  with  the  values 
obtained  for  Alberta  Oil  Incomes  §Z  at  98°C  and  1^8.6 
atmospheres.  The  measured  density  for  the  gas  phase  for 
Alberta  Oil  Incomes  fZ  was  0.1286  g/cc.  Correcting  this 
to  77°  and  167*5  atmospheres  (as  above). gives  a  value  of 
0.1440  g/cc.  The  gas  phase  density  of  31, W.  Hudson* s  Bay 
#6  at  the  same  pressure  and  temperature  is  0.1412  g/cc. 

The  weight  of  pentane  and  heavier  components  per  cc.  of 
gas  phase  is  O.OI58  g  in  U.W.  Hudson’s  Bay  /6  as 
compared  with  0.0207  g  in  Alberta  Oil  Incomes  ^2.  The 
higher  value  in  the  latter  case  can  be  attributed  chiefly 
to  the  much  higher  temperature  and  the  resultant 
increase  in  vapor  pressure  of  the  volatile  liquid 
components.  On  the  basis  of  the  weight  of  pentane  and 
heavier  components  in  the  gas  phase,  it  might  be  expected 
that  the  gas  phase  density  of  Alberta  Oil  Incomes  jj 2 
would  be  greater  than  that  of  N. W.  Hudson’s  Bay  j/6  by 
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0.0049  6/oo.  ie.  a  gas  phase  density  of  0.1461  g/00 
might  be  predicted  for  Alberta  Oil  Incomes  fiz.  This  is 
comparatively  close  to  the  value  of  0.1440  g/co  quoted 
above.  No  account  has  been  taken  of  the  different 
gasroil  ratios  in  comparing  the  two  gas  phase  densities. 
The  lower  gas:oil  ratio  in  the  Alberta  Oil  Incomes  fiZ 
sample  would  be  expected  to  increase  the  gas  phase  density 
but  no  attempt  nas  been  made  to  predict  this  effect  on  a 
quantitative  basis. 

Consideration  of  liquid  phase  densities  shows  that 
increased  pressure  should  cause  a  decrease  in  the  density 
because  of  increased  solubility  of  gaseous  components  at 
the  higher  pressures.  The  increased  transfer  of  volatile 
liquids  from  the  liquid  phase  to  the  gas  phase  with 
increasing  pressure  operates  in  the  reverse  manner  - 
i.e.  to  decrease  the  liquid  phase  density.  This  latter 
effect  is  of  minor  importance  compared  to  the  influence 
of  increasing  solubility  of  gas.  The  decrease  in  liquid 
phase  density  with  increasing  pressure  is  clearly  shown 
in  Runs  1,  2,  J>  and  5»  and  Figure  5. 

An  increase  in  the  gas: oil  ratio  should  cause  an 
increased  liquid  phase  density  due  to  more  complete 
transfer  of  light  liquid  fractions  to  the  gas  phase. 

This,  in  itself,  should  increase  the  liquid  phase  density, 
and  in  addition  there  is  a  decreased  solubility  of  butane 
and  lighter  constituents  in  the  heavier  liquid  residue 
which  is  left  after  transfer  of  the  light  liquids  to  the 
gas  phase.  This  also  would  increase  the  density. 
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Assuming  that  the  approximately  linear  relationship 
between  liquid  phase  density  and  pressure  in  Huns  2, 

3  and  3  continues  to  the  pressure  of  Run  1,  a  liquid 
phase  density  of  0.6377  g/co.  can  be  calculated  for  Run  1. 
The  difference  between  tnis  value  and  the  measured  value 
of  0.6118  can  be  partially  ascribed  to  the  decreased 
gas: oil  ratio  between  Runs  3  and  1. 

The  effect  of  temperature  on  liquid  phase  density 
is  two-fold:  a)  increased  temperature  will  decrease  the 
density  due  to  normal  expansion  of  the  liquid  (2)  increased 
temperature  will  increase  the  density  due  to  increased 
evaporation  of  light  liquid  fractions  into  the  gas  phase 
and  simultaneous  decreasing  solubility  of  gas  in  the  liquid. 

In  Runs  2,  6  there  is  a  small  decrease  in  pressure 
and  a  small  increase  in  gas  roil  ratio,  both  of  which  will 
tend  to  increase  the  liquid  phase  density.  The  total 
effect  will  be  very  small.  The  increase  in  temperature 
will  have  a  two-fold  effect.  The  measured  change  is  a 
decrease  in  liquid  phase  density  which  indicates  that  the 
decrease  in  density  due  to  expansion  with  increasing 
temperature  is  greater  than  the  combined  effects  of  the 
other  three  factors. 

In  Runs  3,  7  there  is  a  slight  decrease  in  pressure 
going  from  77°  to  83°  which  should  increase  the  density; 
simultaneously,  the  small  decrease  in  gas:oil  ratio  should 
decrease  the  density.  Both  effects  will  be  small.  The 
two  opposing  effects  of  temperature  will  likewise  be  felt. 
The  measured  change  in  liquid  phase  density  is  an 


increase.  In  view  of  the  anomalies  previously  discussed 
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in  Run  #7»  it  seems  reasonable  to  assume,  that,  had  a 
state  of  equilibrium  been  attained,  the  liquid  phase 
density  would  have  changed  in  a  manner  analogous  to  the 
ohange  observed  in  the  temperature  interval  70°  -  77°. 

As  explained  previously,  it  was  found  necessary  to 
express  equilibrium  constants  in  terms  of  c  >ncentration 
ratios  rather  than  mole  fraction  ratios.  On  this  basis, 
if  the  system  behaved  ideally,  Henry’s  Law  predicts  a 
constant  K  for  methane  and  for  methane  *  nitrogen. 

The  deviation  from  ideality  which  might  be  expected,  is 
a  decrease  in  K  for  methane  with  increasing  pressure  due 
to  the  more  rapid  increase  in  solubility  of  methane  at  the 
higher  pressures.  The  actual  K’s  for  methane  are  roughly 
constant  for  Runs  1,  2  and  3*  Run  3  is  at  variance  with 
the  others  but  no  explanation  is  suggested.  The  K’s  for 
methane  *■  nitrogen  are  even  more  Constant  except  for 
Run  3  wnich  is  out  of  line.  It  would  seem  that  Henry’s 
Law  does  apply  to  methane  and  methane  *■  nitrogen  in  this 
system,  with  reasonable  accuracy,  within  the  pressure 
range  investigated.  If  the  total  weight  of  methane  from 
both  liquid  and  gas  phase  samples  is  divided  by  the  total 
volume  of  the  two  phase  samples,  a  value  is  obtained  which 
represents  the  approximate  concentration  of  methane  at 
the  critical  pressure.  This  value  is  considerably  lower 
than  the  concentration  of  methane  in  the  gas  phase  and 
higher  than  the  concentration  of  methane  in  the  liquid 
phase  for  all  pressures  investigated.  It  is  thus 
reasonable  to  predict  that  K  for  methane  should  show  a 

decrease  with  increasing  pressure  until  it  reaches  a  value 
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of  1  at  the  critical  pressure  where  the  system  is  single 

phase. 

K  for  ethane  shows  a  continued  increase  with  pressure 
for  all  pressures  investigated.  Here,  again.  Run  {fj,  is 
out  of  line  with  the  rest.  A  calculation  of  the  concentra¬ 
tion  of  ethane  at  the  critical  pressure  shows  that  the 
concentration  of  ethane  should  increase  in  the  gas  phase 
and  decrease  in  the  liquid  phase  with  increasing  pressure. 
Ideally,  K  for  ethane  should  decrease  from  a  large  value 
at  low  pressures  to  a  value  of  1  at  the  critical  pressure. 
In  the  complex  system,  K  for  ethane,  having  apparently 
decreased  to  a  minimum,  is  now  increasing  with  pressure 
and  nas  exceeded  the  value  of  1  at  200  atmospheres.  It 
might  he  predicted  that  a  maximum  value  of  K  will  be 
attained  with  increasing  pressure,  and,  at  the  critical 
pressure,  if  it  could  be  attained,  K  will  have  decreased 
again  to  a  value  of  1. 

Propane  and  butane  both  show  an  initial  decrease  in 
K  with  pressure,  followed  by  an  increase.  This  indicates 
an  increased  transfer  of  these  components  to  the  gas  phase 
at  the  higher  pressures.  This  confirms  the  statement 
made  previously,  that  the  region  of  retrograde  condensation 
is  encountered  in  the  pressure  range  investigated. 

The  general  trend  of  K  for  the  "light  fraction" 
(pentane  *■)  is  to  increase  with  pressure,  but  the 
inaccuracies  in  measuring  this  fraction  are  so  great  tnat 
very  little  significance  can  be  placed  on  the  actual 
numerical  results. 

K  for  methane  **  nitrogen  shows  very  little  change  with 
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temperature.  All  other  K*s  sho^v  a  tandenoy  to  increase 
with  increasing  temperature.  This  oan  be  attributed,  in 
the  case  of  ethane,  to  the  decreased  solubility  of  ethane 
in  the  liquid  phase  at  higher  temperatures.  The  change 
in  K  is  slight.  For  higher  hydrocarbons,  the  increase 
in  K  is  due  to  increased  vaporization  of  these  components 
into  the  gas  phase  at  the  higher  temperatures.  The  effect 
is  very  pronounced  in  tne  case  of  butane. 

The  practical  value  of  these  derived  data  in  their 
application  to  the  Turner  Valley  field  will  now  be 
discussed.  From  the  data,  it  appears  that  the  solubility 
of  butane  and  lighter  in  pentane  and  heavier  components 
at  77°C  and  182.8  atmospheres  is  737*3.5  cu.  ft/bbl.  The 
actual  gas:oil  ratio  of  the  well  operating  at  77°C  and 
173*4  atmospheres  is  approximately  I3OO  cu.  ft/bbl.  This 
indicates  that  II. W.  Hudson^s  Bay  well  was  producing, 
at  the  time  of  sampling,  from  both  a  liquid  and  a  gas 
phase.  This  is  in  direct  contrast  to  Alberta  Oil  Incomes  jfZ 
which  was  considered  to  be  producing  from  a  single  liquid 
phase. 

The  values  listed  under  Go/^p  give  an  indication  of 
the  pentane  ♦*  content  of  the  gas.  This  will  give  the  well 
operator  the  necessary  information  as  to  whether  his  gas 
is  worth  scrubbing  to  recover  these  volatile  fractions. 

Other  figures  in  Table  7  are  of  value  in  applying 
the  so-called  "Brown  plan"  of  operating  wells  in  Turner 
Valley  (2).  This  plan  calls  for  the  removal  of  a 
maximum  of  23  bbls.  of  reservoir  volume  per  acre,  per 
well,  per  day.  From  a  knowledge  of  the  solubility  of  gas 
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in  oil  in  the  reservoir  (Fg/bbl  C^H12f.),  the  gas  produced 
by  the  well  oan  be  divided  into  two  portions:  (1)  gas 
originally  in  solution  in  the  oil  in  the  reservoir, 

(2)  gas  originally  present  as  gas  phase  in  the  reservoir. 
From  the  liquid  phase  density  at  equilibrium,  one  can 
oompute  the  volume  of  liquid  in  the  reservoir  corresponding 
to  1  bbl.  of  oil  produced  at  the  surface.  From  the 
compressibility  of  the  gas  phase,  (Xp/(Vp)g),  one  can 
compute  the  volume  in  the  reservoir  corresponding  to 
1000  ou.  ft.  of  gas  produced  at  the  surface.  With  this 
information  at  hand,  the  operator  needs  only  to  calculate 
the  total  reservoir  volume  corresponding  to  one  barrel  of 
oil  at  the  surface,  at  a  given  gas  roil  ratio  and  from 
this  figure  he  can  readily  compute  his  daily  allowable 
production. 

In  view  of  the  lack  of  calculated  data  on  Alberta 
Oil  Incomes  /2,  it  was  considered  inadvisable  to  attempt 
a  correlation  of  results  with  those  obtained  for  N.W. 
Hudson^  Bay  Such  comparisons  as  have  been  made, 

indicate  a  very  close  similarity  in  the  phase  behavior  of 
the  two  wells.  Since  the  two  wells  are  separated 
geographically  by  almost  the  entire  length  of  the 
producing  horizon,  it  is  apparently  not  an  unjustifiable 
assumption,  to  allot  production  to  all  wells  on  the  same 
basis.  If  the  gas  solubility  in  one  well  had  been  found 
to  be  considerably  greater  than  in  the  other,  then  the 
production  of  the  well  with  the  greater  solubility  might 
reasonably  have  been  based  on  a  larger  reservoir  volume. 

The  fixed  allowable  reservoir  volume  set  by  the  "Brown 
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plan”  seems  to  be  Justified  within  the  limits  of 
engineering  aoouracy. 
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SUMMARY 


The  pressure-volume-temperature  relationships  of 
N..V.  Hudson’s  Bay  #6  natural  gas  have  been  studied 
quantitatively  at  60°C  and  qualitatively  at  2_5°C  and  77°C 
over  the  pressure  range  0  -  27^0  psi.  and  comparison  made 
with  similar  data  for  Royal it e  #27  gas  and  Alberta  Oil 
Incomes  #2  gas.  Residual  volumes  and  fugacities  have 
also  been  calculated  and  comparison  made  with  the  same  two 
gases.  Fugacities  for  N. W.  Hudson’s  Bay  #6  gas  were  also 
compared  with  those  calculated  by  Lewis’  fugacity  rule. 

Phase  equilibrium  studies  in  hydrocarbon  systems 
have  been  continued  and  the  results  of  six  equilibrium 
determinations  have  been  presented.  Tne  measured  and 
derived  data  from  these  det erminat ions  have  been  used 
to  explain  the  behavior  of  the  complex  hydrocarbon 
systems  studied,  and  the  application  of  the  results  to  the 
operation  of  Turner  Valley  has  been  indicated. 
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